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Abstract: 


The  objective  of  this  research  was  to  define  the  conditions  under  which 
pollutants,  in  particular,  those  produced  by  bacteria  such  as  sulfide  end  products  of 
the  SRB,  affect  the  amount  of  hydrogen  absorption  by  iron/steel.  The  research 
focused  on  thiosulfate  and  H2S  as  promoters  of  hydrogen  absorption,  as  well  as 
some  potential  inhibitor  species,  and  included  both  modeling  and  hydrogen 
permeation  experiments.  For  improved  understanding  of  how  these  species  affect 
the  absorption  of  hydrogen  into  steel,  the  lyer-Pickering-Zamenzadeh  (IPZ) 
analysis  method  was  first  generalized  to  include  Frumkin  adsorption  of  the 
hydrogen  on  the  steel’s  surface  and  then  used  to  analyze  the  measured  steady  state 
hydrogen  permeation  data.  The  results  are  described  in  the  six  publications  to  date 
and  one  review  paper  that  is  in  preparation. 


Overview 


The  absorption  of  electrolytic  hydrogen  into  the  steels  used  in  the  marine 
industry  leads  to  their  degradation  due  to  hydrogen  embrittlement.  It  is  such  a 
serious  problem  that  the  hydrogen  embrittlement  of  steels  has  been  a  major  focus 
of  researchers  for  several  decades.  It  is  generally  accepted  that  once  hydrogen  is 
absorbed,  the  steel  will  be  susceptable  to  one  form  or  another  of  hydrogen 
degradation.  The  best  known  examples  are  those  encountered  in  environments 
where  hydrogen  sulfide  is  present  such  as  in  the  presence  of  micro-organisms,  e.g., 
sulfate  reducing  bacteria  (SRB)  which  produce  hydrogen  sulfide  in  the  cycle  of 
their  metabolism. 

Hydrogen  absorption  into  metals  occurs  as  a  side  reaction  of  the  hydrogen 
evolution  reaction  (HER)  on  metals.  The  HER  occurs  as  a  viable  cathodic 
reaction  during  corrosion,  cathodic  protection,  pickling  and  electroplating  of 
metals.  It  is  because  of  the  multi-step  nature  of  the  HER  that  hydrogen  absorption 
occurs. 

The  detailed  results  of  this  investigation  are  described  in  the  following 
published  research  papers  and  one  review  paper  that  is  in  preparation. 

“Effect  of  Thiosulfate  and  Sulfite  on  the  Permeation  Rate  of  Hydrogen  Through  Iron”, 
Corrosion,  57, 428  (2001):  M.H.  Abd  Elhamid,  B.G.  Ateya,  K.G.  Weil,  and  H.W. 
Pickering 

“An  Analysis  Procedure  for  Hydrogen  Absorption  under  Frumkin  Adsorption 
Conditions”,  J.  Electrochemical  Soc.,  148,  E248  (2001):  F.M.  Al-Faqeer  and  H.W. 
Pickering 

“Determination  of  the  Rate  Constants  of  Hydrogen  Absorption  into  Metals”,  J. 
Electrochemical  Soc.,  147,  2959  (2000):  M.H.  Abd  Elhamid,  B.G.  Ateya  and  H.W. 
Pickering 

“Calculation  of  the  Hydrogen  Surface  Coverage  and  Rate  Constants  of  the  Hydrogen 
Evolution  Reaction  from  Polarization  Data”,  J.  Electrochemical  Soc.,  147,  2148  (2000): 
M.H.  Abd  Elhamid,  B.G.  Ateya,  K.G.  Weil  and  H.W.  Pickering 

“Effect  of  Benzotriazole  on  the  Hydrogen  Absorption  by  Iron”,  J.  Electrochemical  Soc. 
144,  L58  (1997):  M.H.  Abd  Elhamid,  B.G.  Ateya  and  H.W.  Pickering 

‘The  Effect  of  Iodide  Ions  on  the  Kinetics  of  Hydrogen  Absorption  by  Iron”,  J. 
Electrochemical  Soc.,  147,  2258  (2000):  M.H.  Abd  Elhamid,  B.G.  Ateya  and  H.W. 
Pickering 

“Kinetics  of  Electrochemical  Hydrogen  Degradation  in  Metals”  in  preparation:  F.M.  Al- 
Faqeer  and  H.W.  Pickering 


KINETICS  OF  ELECTROCHEMICAL  HYDROGEN  DEGRADATION  IN  METALS 


F.  M.  Al-Faqeer  and  H.  W.  Pickering 

Department  of  Materials  Science  and  Engineering 
The  Pennsylvania  State  University 
University  Park,  PA  16802 


Abstract 

This  paper  will  review  the  recent  developments  regarding  the  kinetics  of  hydrogen  entry  into 
metals.  Mechanistic  analyses  have  been  developed  which  evaluate  the  kinetics  of  the  hydrogen 
evolution  and  absorption  reactions  on  metals.  These  analyses  can  determine  the  individual  rate 
constants  of  the  hydrogen  evolution  reaction,  k\  and  k2,  (from  the  steady  state  permeation  data  or 
the  polarization  curve)  and  the  kinetic-diffusion  rate  constant,  k,  the  hydrogen  surface  coverage, 
0H,  and  the  hydrogen  concentration  in  the  charging  side  of  the  membrane,  C° ,  (from  permeation 
data).  The  original  IPZ  analysis  has  been  used  extensively  in  the  literature.  An  example  is 
given  in  this  review  illustrating  the  application  of  the  newer  generalized  IPZ  analysis  for 
characterizing  the  effect  of  hydrogen  sulfide  on  the  hydrogen  evolution  and  absorption 
reactions. 
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Introduction 

The  hydrogen  absorption  and  evolution  reactions  on  metals  have  been  studied  extensively  in  the 
literature  after  the  development  of  the  hydrogen  electrochemical  permeation  cell*"’®.  Progress  in 
combining  the  hydrogen  absorption  reaction  with  the  hydrogen  evolution  reaction  was  made  in 
different  studies  ’  ’  °  Different  analyses  have  been  used  in  attempts  to  evaluate  the  kinetics 
of  the  electrochemically  charged  hydrogen  permeation  through  metals.  This  paper  will  focus  on 
the  IPZ  analyses,  the  first  of  which  was  developed  by  Iyer  et.  al.,  in  1989  to  evaluate  the 
hydrogen  permeation  data  and  to  describe  the  kinetics  of  the  hydrogen  evolution  and  absorption 
reactions  . 


The  Original  IPZ  Analysis 

The  original  IPZ  analysis  developed  a  set  of  equations  which  relate  the  charging,  recombination 
and  steady  state  hydrogen  permeation  current  densities  and  the  charging  electrode  potential. 
The  analysis  has  been  used  widely  to  analyze  the  hydrogen  permeation  data  on  iron  in  the 
presence  of  different  solution  additives,  inhibitors  and  promoters"''^’’^.  The  analysis  has  also 
been  used  to  analyze  the  permeation  data  for  different  metals' 

The  hydrogen  evolution  reaction  (HER)  is  not  a  single  step  reaction.  The  reaction  proceeds  in 
two  main  steps,  one  of  which  is  the  discharge  of  protons  on  the  metal  surface  to  form  adsorbed 
hydrogen  atoms.  The  discharge  step  is  represented  by  reaction  (1). 


+  e  +  M->  M-Hads 


(  1  ) 


In  a  corrosive  environment,  the  majority  of  the  adsorbed  atomic  hydrogen  (M-Hads)  combines  to 
form  molecular  hydrogen  which  bubbles  off  the  surface.  This  can  occur  by  either  of  two 
reactions,  but  the  IPZ  analysis  only  applies  when  it  occurs  by  chemical  recombination  as  shown 
in  reaction  (2). 


M-Hads  +  M-Hads  -^2M  +  H2 


(  2  ) 


A  side  reaction  of  the  HER  is  the  hydrogen  absorption  reaction  (HAR)  which  takes  place 
according  to  reaction  (3). 


M-Hads  -  M  +Habs 


(  3  ) 


where  Habs  refers  to  a  hydrogen  atom  absorbed  in  the  metal. 

The  IPZ  analysis  assumes  that  the  HER  proceeds  according  to  a  coupled  discharge- 
recombination  (Volmer-Tafel)  mechanism.  In  the  original  IPZ  analysis  Langmuir  conditions 
were  assumed  for  the  hydrogen  surface  coverage.  Thus,  the  analysis  ignores  the  adsorption  of 
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other  species  on  the  surface.  The  analysis  also  assumes  that  the  charging  current  is  the  sum  of 
the  recombination  and  steady  state  permeation  current  densities  (ic  =  ir  +  i«.)'  '• 

The  original  IPZ  analysis  expresses  the  current  density  of  the  first  step  of  the  HER  (reaction  1) 
using  the  Langmuir  adsorption  isotherm  as  follows'* 


ic  =  F/ciCh+  (  1-6h)  exp(  -aar]) 


(  4  ) 


where  4  is  the  cathodic  current  density,  Ch+  is  the  hydrogen  ion  concentration,  a  =  F/RT,  a  is 
the  transfer  coefficient  and  r|  is  the  hydrogen  overpotential. 

Equation  (4)  was  rearranged  in  the  analysis  to 


4exp(  aar|)  =  4'  (  1-0h  ) 

(  5  ) 

4'  =  F^iCH+=4/(l-e'H) 

(  6  ) 

where  4exp(aaTi)  is  referred  to  as  the  charging  fimction,  4  is  the  exchange  current  density  of  the 
hydrogen  evolution  reaction  and  6^  is  the  equilibrium  hydrogen  surface  coverage. 

The  rate  of  the  recombination  step  of  the  hydrogen  evolution  reaction  (reaction  2)  was  given  by 
the  analysis  as"’*^ 


4  =  Fk2Q\ 


(  7  ) 


The  flux  of  the  hydrogen  atoms  passing  through  the  membrane  using  Pick's  first  law  is  shown  in 
equation  (8).^ 


J  =  -D 


dx 


(  8  ) 


Using  equations  (9)  and  (10)  as  the  boundary  conditions, 


C(x  =  0)=C“ 


(  9  ) 


C(x  =  Z)=  =  0 


(  10  ) 
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the  steady  state  hydrogen  permeation  flux  is’*^ 


C° 

J  =  D^ 
L 


(  n  ) 


where  D  is  the  hydrogen  diffusion  coefficient  within  the  metal  and  L  is  the  thickness  of  the 
membrane. 

From  the  flux  in  equation  (11),  the  steady  state  hydrogen  permeation  current  density  is  given 

by  >’2 


C" 

”  L 


(  12  ) 


Equation  (12)  can  be  used  to  estimate  the  absorbed  hydrogen  concentration  at  the  charging  side 
of  the  membrane,  ,  if  tire  hydrogen  diffusivity  in  the  membrane  is  known.  The  diffusion 
coefficient  of  hydrogen  through  the  membrane  can  be  found  using  the  permeation  transient  firom 
the  breakthrough  time,  tb.  The  breakthrough  time  is  the  time  it  takes  the  hydrogen  atom  to 
diffuse  through  the  membrane  from  the  charging  to  the  exit  side.  It  can  be  calculated  using 
equation  (13).'* 


h 


15.3*X> 


(  13  ) 


Kim  and  Wilde*,  showed  that  the  steady  state  hydrogen  permeation  current,  /«>,  is  given  by 
equation  (14).  The  IPZ  analysis  uses  this  relation  in  the  development  of  its  main  relations". 


P^des^H 


(  14  ) 


Substituting  from  equation  (12)  into  equation  (14)  gives 


F^des 


LL 

_ CO 

VD 


Rearranging  equation  (15)" 


(  15  ) 
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L=F 


''abs 


1+^ 


D 


(  16  ) 


The  kinetic-diffusion  constant,  k,  is  expressed  in  terms  of  the  adsorption  and  absorption  rate 
constants  as  shown  in  equation  (17)"’^^’ 


k  = 


''abs 


l  +  ^des  — 

des  ^ 


(  17  ) 


Equation  (16)  can  be  simplified  using  the  kinetic-diffusion  constant,  k,  to  equation  (18). 


=  F  ^  9h 


(  18  ) 


Equation  (18)  relates  the  steady  state  hydrogen  permeation  current  density  to  the  hydrogen 
surface  coverage.  Equations  (17)  and  (18)  are  two  major  relations  in  the  IPZ  analysis. 

Substituting  the  hydrogen  surface  coverage  from  equation  (7)  into  equation  (18)  and  rearranging 


L  =  k 


(  19  ) 


Equation  (19)  shows  the  relation  between  the  steady  state  hydrogen  permeation  current  density, 
ice,  and  the  square  root  of  the  hydrogen  recombination  current  density,  ^ .  This  relation  should 


be  a  straight  line  passing  through  the  origin  with  a  slope  of  k  I —  for  the  application  of  the  IPZ 


analysis.  This  limitation  on  the  application  of  the  original  IPZ  analysis  will  be  discussed  in 
more  detail  below. 


Similarly,  substituting  the  hydrogen  coverage,  Gj^,  from  equation  (18)  into  equation  (5)  and 
rearranging  gives". 


/cexp(  aari)  =  zV  (1  ~  L) 


(  20  ) 
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Equation  (20)  shows  the  relation  between  the  charging  function,  Zcexp(aari)  and  the  steady  state 
hydrogen  permeation  current  density,  This  relation  should  be  a  straight  line  with  a  slope  of 

and  an  intercept  of  /' . 

The  original  IPZ  analysis  developed  the  major  relations,  fl7),  (18),  (19)  and  (20),  between  the 
various  variables  in  the  hydrogen  permeation  system’ The  kinetic-diffusion  constant, 
k,  and  4'  can  be  found  from  the  slope  and  intercept  of  equation  (20).  The  recombination  rate 
constant,  ki,  can  be  evaluated  using  the  estimated  Wnetic-diffusion  constant,  k,  and  the  slope  of 
equation  (19).  The  value  of  the  exchange  current  density,  io,  can  be  estimated  using  the  value  of 
/o'.  The  discharge  rate  constant,  k\,  can  be  found  from  the  exchange  current  density  using  the 
formula  /o'  =  FA:iCh.  The  hydrogen  surface  coverage,  0h,  can  be  evaluated  using  the  kinetic- 
diffusion  constant,  k,  and  equation  (18).  Therefore,  these  kinetic  and  thermodynamic 
parameters  of  both  the  hydrogen  evolution  and  absorption  reactions  can  be  determined  using  the 
IPZ  analysis  when  the  experimental  data  satisfies  the  requirements  of  equations  (19)  and  (20). 


Membrane  Thickness  Effect 

The  IPZ  analysis  can  be  used  to  evaluate  the  absorption  and  desorption  rate  constants ,  itabs  and 
Ades-  Rearranging  equation  (17)  can  give  a  relation  between  the  kinetic-diffusion  constant,  k,  and 
the  membrane  thickness,  L. 


1  = 


(  21  ) 


Permeation  experiments  need  to  be  done  at  different  metal  membrane  thicknesses,  L,  in  order  to 
use  equation  (21).  Once  the  kinetic-diffusion  constants  are  obtained  at  different  membrane 
thicknesses,  equation  (21)  can  be  used  to  obtain  Aabs  and  Araes-  Equation  (21)  shows  a  straight  line 
between  the  reciprocal  of  the  kinetic-diffusion  constant,  k,  and  the  membrane  thickness,  L.  The 
slope  of  this  relation  will  give  the  value  of  the  desorption  constant,  ^des-  The  intercept  will  give 
the  absorption  rate  constaint,  Arabs- 


Calculation  of  9h  and  rate  constants  of  the  HER  from  Polarization  Data 

Elhamid  et  al.  (21),  showed  that  the  polarization  curve  of  the  HER  can  be  analyzed  to  calculate 
the  hydrogen  surface  coverage  and  the  rate  constants  of  the  hydrogen  discharge  and 
recombination  reactions  for  metals  with  very  low  hydrogen  permeability.  They  assumed  in  their 
study  that  the  permeation  current  density  is  very  low  for  some  metals  such  as  copper  and 
therefore,  the  charging  current  density,  ic,  will  be  equal  to  the  recombination  current  density,  /,. 
They  developed  and  used  Equation  22  to  calculate  the  exchange  current  density,  /q,  and  Aa. 
Then,  ki  was  calculated  from  the  exchange  current  density  using  io  =  FAiCh+.  The  hydrogen 


surface  coverage  was  calculated  from  Equation  7.  More  details  about  this  kind  of  calculation 
can  be  found  in  Ref  21. 
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/cexp(  aari)  =  4'  (1 


(  22  ) 


Frumkin  Modification  of  the  Original  IPZ  Analysis 

The  Frumkin  adsorption  isotherm  assumes  that  the  free  energy  of  adsorption  of  a  species 
decreases  with  coverage  according  to  the  following  equation^^'^^. 


AGe“  =AG°+/R7’9 


(  23  ) 


where  /  is  a  dimensionless  factor  that  describes  the  deviation  from  the  ideal  Langmuir  behavior. 
The  /  RT  quantity  is  sometimes  called  the  rate  of  change  of  the  standard  free  energy  of 

adsorption,  AG^  and  AGq  are  the  standard  free  energy  of  adsorption  at  a  certain  coverage,  9, 
and  at  zero  coverage,  0  =  0,  respectively,  R  is  the  gas  constant  and  T  is  the  absolute  temperature. 

The  Langmuir  adsorption  isotherm  assumes  that  the  standard  free  energy  of  adsorption  is 
independent  of  coverage  (/  =  0).  This  isotherm  is  usually  applicable  at  low  values  of  coverage 
(0  <  0.1)  and  at  coverages  close  to  unity  (0  >  0.9)^'^’^^  In  some  instances  it  has  been  found  to  be 
applicable  also  for  intermediate  coverages  of  hydrogen  on  iron  surfaces. 

Writing  the  rate  equations  in  view  of  the  Frumkin  adsorption  conditions,  one  has  for  the  rate  of 
proton  discharge  and  hydrogen  evolution. 


4  =  ^'o  ( 1  ~  0  H )  exp(-aQr  ?j)exp(-a  /  0  „  )  (  24  ) 

i'r  =  FA^jGjj  exp(2Qr  0jj)  (  25  ) 


Inserting  0h  from  equation  (1 8)  into  equation  (25) 


(  26  ) 


Taking  the  square  root  of  equation  (26) 
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V7  =  ^ ‘expC,./^) 


(  27  ) 


Dividing  both  sides  of  equation  (27)  by  z’oo  and  taking  the  natural  logarithm  (In)^^*^^ 


ln(^)  =  ln(M  *  7  )  +  ^;, 

V  F  k  Fk  " 


(  28  ) 


This  relation,  which  was  also  derived  by  Iyer,  et.  al.'^’^*,  in  a  similar  form,  but  with  different 
slope  and  intercept,  can  be  used  for  determining  the  rate  constants  when  the  value  off  is  known 
or  can  otherwise  be  estimated.  Recently,  Al-Faqeer  and  Pickering  (29)  arrived  at  another 
relationship  which,  when  applied  to  the  steady  state  permeation  data,  can  yield  the  rate  constants 

wiAout  a  prior  knowledge  of  the  value  off.  In  fact,  the  value  of  /  is  a  result  of  the  analysis 
This  derivation  follows: 

Rearranging  equation  (24)  and  solving  for  exp(a/9H) 


exp(a  / e„ )  =  ^h) 

4exp(aa;7) 


(  29  ) 


Taking  the  square  root  of  equation  (25) 


7^  =  7^  *  0fj  *exp(a/0j,) 


(  30  ) 


Substituting  exp(<2r  /0fj)  from  equation  (28)  into  equation  (30) 


♦Ao-^h) 

4exp(a  a  tj) 


(  31  ) 


Rearranging  and  inserting  for  0h  from  equation  (18)  into  equation  (3 1) 


7^  4exp(a  a  7)  =  7f^  *  ^  (1  _  ^) 

tk  Fk 


(  32  ) 
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Dividing  both  sides  of  equation  (32)  by 


4exp(aa?7)  = 


(  33  ) 


Equations  (28)  and  (33)  can  be  simplified  to  equations  (19)  and  (20)  if  tlie  Langmuir  adsorption 
conditions  are  assumed  instead  of  the  Frumkin  adsorption  conditions.  In  other  words,  the 
generalized  IPZ  analysis^^  reduces  to  the  original  Langmuir  IPZ  analysis^’  if  the  standard  free 
energy  of  adsorption  is  assumed  to  be  independent  of  the  hydrogen  coverage,  i.e.,  /  =  0.  Thus, 
the  original  IPZ  analysis  is  a  special  case  of  the  generalized  IPZ  analysis,  in  the  same  way  the 
Langmuir  isotherm  is  a  special  case  of  the  Frumkin  isotherm. 


Plots  of  both  In(^^p-)  and  i^exp(aa  tf)  vs.  /«  should  give  straight  lines  according  to 


equations  (28)  and  (33),  respectively. 


The  procedure  of  calculating  the  unknowns  in  the  system  is^^: 

>  From  the  slope  and  intercept  of  equation  (33),  k  can  be  calculated. 

>  From  the  slope  of  equation  (28)  and  knowing  k  from  step  1 ,  /  can  be  calculated,  after  using 
the  charging  (Tafel)  data  to  obtain  the  value  of  or. 

>  From  the  intercept  of  equation  (28)  and  knowing  k  from  step  1 ,  ^2  can  be  calculated. 

^  From  the  intercept  of  equation  (33)  and  knowing  k  from  step  1  and  kz  from  step  3,  can  be 
calculated. 

^  From  equation  (18)  and  knowing  k  from  step  1,  the  hydrogen  surface  coverage  (0h)  can  be 
calculated. 


Application  of  the  Generalized  IPZ  analysis 

The  original  IPZ  analysis  has  been  used  extensively  in  the  literature.  Details  about  using  the 
original  IPZ  analysis  to  evaluate  the  experimental  data  can  be  found  elsewhere'  '>19,20,28^  what 
follows,  we  give  an  example  where  the  Frumkin  modification  of  the  IPZ  analysis  was  used  to 
evaluate  the  experimental  data  on  the  effect  of  hydrogen  sulfide  on  the  hydrogen  evolution  and 
absorption  reactions  on  an  iron  membrane.  The  electrochemical  hydrogen  permeation  technique 
was  used  to  carry  out  the  experimental  part  of  the  research. 

The  effect  of  hydrogen  sulfide  concentrations  on  the  relation  between  the  measured  steady  state 
hydrogen  permeation  current  density,  /«,,  and  the  cathodic  potential,  E,  in  the  charging 
compartment  of  the  cell  is  presented  in  Figure  1.  The  figure  shows  that  hydrogen  sulfide 
strongly  enhances  hydrogen  absorption  into  iron  since  the  steady  state  hydrogen  permeation 
current  density  increases  by  more  than  an  order  of  magnitude  with  increasing  H2S  concentration 
at  the  same  cathodic  potential. 
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Figure  2  shows  the  effect  of  hydrogen  sulfide  concentrations  on  the  relation  between  the  rate  of 
the  hydrogen  charging  reaction,  z'c,  and  the  cathodic  potential,  E.  The  figure  reveals  that 
hydrogen  sulfide  enhances  the  HER  since  the  hydrogen  overpotential  is  less  {E  is  more  positive) 
at  the  same  value  of  the  applied  cathodic  charging  current. 


Figure  1:  The  relation  between  /«,  and  E  obtained  at  different  hydrogen  sulfide  concentrations  at 
pH  2 


Figure  2:  The  relation  between  4  and  E  obtained  at  different  hydrogen  sulfide  concentrations  at 
pH  2 


Figure  3  shows  the  effect  of  hydrogen  sulfide  on  the  hydrogen  concentration  inside  the 
membrane.  This  figure  reveals  that  the  hydrogen  concentration  inside  the  membrane  increases 
as  the  hydrogen  sulfide  concentration  increases  in  the  charging  solution. 
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Figure  3:  The  relation  between  C°  and  E  obtained  at  different  hydrogen  sulfide  concentrations 
at  pH  2 


Figure  4  shows  the  relation  between  the  measured  steady  state  permeation  current  density, 
and  the  square  root  of  the  recombination  current  density,  ,  obtained  at  different  hydrogen 

sulfide  concentrations.  This  Figure  shows  straight  lines  but  only  the  blank  solution  line  passes 
by  the  origin  in  accord  with  equation  (19).  In  the  presence  of  hydrogen  sulfide,  the  lines  do  not 
pass  through  the  origin  and  have  a  finite  intercept  with  the  y-axis.  This  indicates  that  the 
original  IPZ  analysis  cannot  be  used  to  evaluate  the  data  in  the  presence  of  hydrogen  sulfide. 


Figure  4:  The  relation  between  /«,  and  ^  obtained  at  different  hydrogen  sulfide 
concentrations  at  pH  2. 
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Figure  5  shows  the  relation  between  ]n(-^-^)  and  the  measured  steady  state  hydrogen 


permeation  current  density,  /«,.  The  figure  shows  straight  lines  in  accord  with  equation  (27). 
Figure  6  shows  linear  relations  between  /  exp(ao:;7)  and  the  measured  steady  state 


hydrogen  permeation  current  density,  i„,  in  accord  with  equation  (32).  This  indicates  the 
generalized  IPZ  analysis  can  be  used  to  analyze  the  data. 


Figure  5  The  relation  between 
concentrations  at  pH  2 


ln(^) 


and  im  obtained  at  different  hydrogen  sulfide 


The  slopes  of  the  lines  in  Figure  6  were  used  along  with  the  intercepts  for  the  values  of  the 
kinetic-diffusion  constant,  k,  according  to  Equation  (32).  The  values  of  k  at  different  hydrogen 
sulfide  concentrations  were  used  to  estimate  the  hydrogen  recombination  rate  constant,  k2,  and 
/,  from  the  intercepts  and  slopes  of  the  lines  in  Figure  6  respectively,  in  accord  with  Equation 
(27).  Values  of  k  and  k2  were  used  to  calculate  the  hydrogen  exchange  current  density  at 

different  hydrogen  sulfide  concentrations,  —  =  i^  when  9^«1,  from  the  intercepts  of 

(1“9h) 

Figure  6,  in  accord  with  equation  (32).  The  calculated  values  of  the  exchange  current  density 
were  used  to  estimate  the  discharge  rate  constant,  k\,  at  different  hydrogen  sulfide  concentration 
using  the  relation  io=Vk\Cn+. 
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Figure  6: 


The  relation  between 


/(.expCaari)  and  obtained  at  different  hydrogen  sulfide 


concentrations  at  pH  2 


Table  I  shows  values  for  i^,  ku  ki,  k  and  /  at  different  hydrogen  sulfide  concentrations.  These 
results  show  that  hydrogen  sulfide  increases  the  proton  discharge  rate  constant,  kx,  and  decreases 
the  recombination  rate  constant  of  the  HER,  kj.  This  could  lead  to  an  increase  in  the  hydrogen 
surface  coverage,  0h.  The  values  of  k  were  used  to  estimate  the  hydrogen  surface  coverage,  0h, 
using  Equation  (18)  at  different  potentials  and  hydrogen  sulfide  concentrations.  Figure  7  shows 
the  relations  between  the  hydrogen  surface  coverage,  0h,  and  cathodic  potential,  E.  It  reveals 
that  the  hydrogen  coverage  increases  with  the  increase  in  both  hydrogen  sulfide  concentration 
and  the  potential  in  the  cathodic  direction. 


E,  V(SCE) 


Figure  7:  The  relations  between  0h  and  E  obtained  at  different  hydrogen  sulfide  concentrations 
at  pH  2 


13 


Table  I  also  shows  that  hydrogen  sulfide  affects  f.  The  value  of  /  increased  from  almost  zero 
(0.49)  in  the  absence  of  hydrogen  sulfide  to  nearly  8.7  in  the  presence  of  hydrogen  sulfide.  This 
indicates  that  the  standard  free  energy  of  adsorption  is  coverage  dependent,  and  that  the 
Frumkin  adsorption  isotherm  was  correctly  applied. 

The  increase  in  0h  and  /,  dominate  over  the  decrease  in  the  hydrogen  recombination  rate 
constant  of  the  HER,  k2,  which  leads  to  an  increase  in  the  hydrogen  evolution  reaction  rate. 
However,  the  increase  in  the  hydrogen  absorption  reaction  rate  is  attributed  to  increases  in  the 
hydrogen  surface  coverage,  Gh,  and  the  kinetic-diffusion  constant,  k. 


Table  I:  Values  of  k,  kj,  io,  h  and  f  obtained  on  iron  membrane  at  different  hydrogen  sulfide 

concentrations  at  pH  2. 


HzS 

Content 

k 

mol  cm'^  s‘^ 

k2 

mol  cm'^  s'* 

^*0 

p  A  cm'^ 

h 

cm  s'* 

/ 

Blank 

7.02X  10‘“ 

4.96  X  lO-'^ 

0.421 

4.36x10'* 

0.49 

50  pM 
(1.7  ppm) 

6.61  X  lO'*'^ 

2.82  X  10'" 

2.78 

2.89  X  10'" 

8.94 

500  pM 
(17  ppm) 

9.62  X  lO'*^ 

2.66  X 

9.22 

9.56  X  10'" 

8.87 

5mM 
(170  ppm) 

1.43  X  10'^ 

9.24  X  10'*" 

26.1 

2.70  X  10'^ 

8.26 

Conclusion 

The  kinetics  of  hydrogen  entry  in  metals  has  been  studied  using  different  analyses.  The  original 
IPZ  analysis  has  been  applied  to  steady  state  permeation  data  to  yeild  the  individual  rate 
constants  of  the  hydrogen  evolution  reaction,  k\  and  k2,  the  kinetic-diffusion  constant,  k,  the 
hydrogen  surface  coverage,  0h,  and  the  hydrogen  concentration  in  the  charging  side,  C°.  It  was 
also  shown  how  0h,  k\  and  k2  could  be  determined  from  the  measured  HER  polarization  curve. 
The  analysis  was  modified  assuming  the  validity  of  the  Frumkin  adsorption  isotherm  where  a 
general  analysis  was  developed.  The  experimental  data  in  the  presence  of  hydrogen  sulfide  fits 
the  relations  derived  in  the  case  of  tlie  generalized  IPZ  analysis.  However,  there  is  no 
independent  evidence  that  the  Frumkin  approach  is  a  better  description  than  the  Langmuir 
approach  for  hydrogen  adsorption  under  these  conditions.  Furthermore,  the  adsorption  of  the 
sulfur  species  was  ignored;  therefore,  any  effect  that  they  may  have  had  on  the  available 
hydrogen  adsorption  sites  on  the  surface  was  ignored  too,  for  the  sake  of  an  approximation. 
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List  of  Symbols 


a  a  constant,  F/RT,  V' 

Ch^  hydrogen  ion  concentration  in  the  electrolyte,  mol  L*' 

C°  hydrogen  concentration  in  the  metal  at  the  charging  surface,  ppm 

D  hydrogen  diffusion  coefficient,  cm^  s'' 

F  Faraday  constant,  96484.6  C  mof' 

/  a  factor  that  describes  the  deviation  from  the  ideal  Langmuir  behavior, 

dimensionless 

AGg°  standard  free  energy  of  adsorption  at  certain  coverage,  J  mof' 

AGq  standard  free  energy  of  adsorption  at  zero  coverage,  J  mof' 

ic  charging  current  density,  A  cm'^ 

Zoo  steady  state  hydrogen  permeation  current  density,  A  cm’ 

z’o  exchange  current  density,  A  cm'^ 

z’o'  Acm'^ 

z'r  hydrogen  recombination  current  density,  A  cm'^ 

Arabs  absorption  rate  constant,  mol  (cm^  s)'' 

Ardes  desorption  rate  constant,  cm  s"' 

k  kinetic-diffusion  constant,  mol  cra'^  s’' 

Ati  discharge  rate  constant,  cm  s*' 

k2  recombination  rate  constant,  mol  (cm^  s)’' 

L  membrane  thickness,  cm 

R  gas  constant,  8.3 14  J  (mol  K)’' 

T  absolute  temperature,  K 

Greek  Symbols 

a  transfer  coefficient,  dimensionless 

T]  overpotential  of  the  HER,  V 

9ii  hydrogen  surface  coverage,  dimensionless 

hydrogen  surface  coverage  at  equilibrium,  dimensionless 
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ABSTRACT 

The  results  presented  show  that  benzotriazole  (BTA)  Inhibits  the  absorption  of  hydrogen  into  Iron  which  Is  cathodically  polarized  in  an 
acid  (pH  3  1 .7)  sulfate  solution.  Analysis  of  the  results  indicates  that  BTA  Inhibits  also  the  hydrogen  evolution  reaction,  although  It  does  ; 
not  change  Its  mechanism,  which  Is  shown  to  be  proton  discharge-Tafel  recombination.  However,  BTA  shifts  the  position  of  the  equllibrl-  j 
urn,  H,^  H,^„  toward  the  left  side  and  hence  leads  to  a  decrease  In  the  concentration  of  within  the  lattice.  The  extent  of  this  | 

decrease  depends  on  the  BTA  concentration,  e.g..  about  fivefold  at  1  mM  and  tenfold  at  1 0  BTA.  This  effect  is  opposite  to  that  report-  1 

ed  for  Inhibitors  such  as  thiourea  which  inhibits  the  hydrogen  evolution  reaction  and  yet  promotes  hydrogen  permeation. 


Introduction 

The  absorption  of  hydrogen  In  metals  and  alloys  Is  a  precursor 
to  hydrogen  embrittlement,  the  deterioration  of  the  materials 
mechanical  properties,  and  failure  by  crack  propagation.**  The 
problem  is  hardly  avoidable  since  hydrogen  enters  metals  and 
alloys  from  many  sources,  eg,,  corrosion,  acid  pickling,  electro¬ 
plating,  welding,  cathodic  protection.  In  aqueous  environments, 
the  introduction  of  hydrogen  within  the  metal  or  the  alloy  occurs 
during  the  hydrogen  evolution  reaction  (HER).’*  There  are  two 
generally  accepted  mechanisms  for  this  reaction’*^  the  dis- 
chafge-Tafel  recombination 


H'  -r  e  -M-H,^  [1] 

M-H„^  +  -  Hj  |2] 

and  the  discharge-electrochemical  desorption  which  begins  with 
reaction  1  followed  by 


H‘  -I-  e  [3] 

where  M-H,^  refers  to  a  nascent  hydrogen  atom  adsorbed  on  the 
metal  surface.  The  major  part  of  this  adsorbed  hydrogen  reacts  to 
give  hydrogen  molecules  through  reactions  2  or  3.  A  fraction  of  this 
adsorbed  hydrogen  is  adsorbed  Into  the  metallic  lattice, 

M-H^  ^  M(H,b,)  [4] 

This  adsorbed  hydrogen  diffuses  within  the  metal  and  eventually 
causes  embrittlement.  Many  theories  have  been  proposed  to 
explain  the  pheniomenon,  e.g.,  the  pressure  theory,  the  adsorption 
theory,  the  decoliesion  theory,  and  the  hydride  theory.*  * 

Different  approaches  have  been  proposed  to  mitigate  tlie  prob¬ 
lem,  either  by  the  use  of  inhibitors,®  by  coating  the  surface.’* or 
alloying  with  elements  which  have  lower  diffusivilies  for  hydro¬ 
gen.*®*®  Among  the  inhibitors,  BTA  has  been  tested  on  iron*  and 
steel*®  in  strongly  acidic  media.  None  of  these  referenced  works 
dealt  with  the  mechanism  Involved  in  the  inhibition  of  the  hydrogen 
absorption  process.  Further,  both  sets  of  results  produced  different 
conclusions.  For  Instance,  Dull  and  Nobe*  found  tfiat  the  efficiency 
of  BTA  as  an  inhibitor  for  hydrogen  absorption  Increases  signifi¬ 
cantly  with  the  rate  of  hydrogen  evolution,  whereas  Subramanyan 
et  a/.’°  found  that  the  efficiency  to  be  independent  of  the  extent  of 
cathodic  polarization  and  hence  of  the  rate  of  hydrogen  evolution. 
These  measurements  were  obtained  in  unbuffered  and  strongly 
acidic  media  where  BTA  is  known  to  be  less  efficient  as  a  general 
corrosion  inhibitor.** 

Our  objective  here  is  to  explore  the  effect  of  BTA  on  tfie  kinetics 
of  the  hydrogen  absorption  within  iron.  The  experimental  results 
were  obtained  in  a  buffered  acidic  medium  of  a  pH  at  which  BTA 
has  a  high  efficiency  as  a  general  corrosion  inhibitor.  BTA  is  wide¬ 
ly  used  as  a  general  corrosion  inhibitor  for  copper  and  its  alloys** 
and  to  lesser  extent  for  iron  base  alloys.*'**** 


•  Electrochemical  Society  Student  Member. 

**  Electrochemica!  Society  Active  Member. 

**  Elecirocfiemical  Society  Fellow. 


Experimental 

A  hydrogen  permeation  celP’®  was  used  to  measure  the  per¬ 
meability  of  hydrogen  through  iron  membranes.  The  cell  consists 
of  two  Identical  parts  separated  by  the  metal  membrane,  Fig.  1. 
One  side  and  Its  associated  circuitry  were  used  to  generate  the 
hydrogen  charging  current  and  to  measure  the  potential  of  the 
charging  side  of  the  membrane  while  the  other  side  was  used  to 
monitor  the  diffusion  flux  of  hydrogen  through  the  membrane  by 
measuring  the  hydrogen  oxidation  current  at  the  exit  surface  of  the  ! 
membrane.  The  charging  side  had  a  saturated  calomel  reference 
electrode  whereas  the  exit  side  had  a  Hg/HgO  reference  electrode. 
Both  sides  had  Pt  counlerelectrodes,  which  were  separated  from 
the  solution  using  glass  frits.  The  exposed  surface  area  of  the 
sample  was  0.8  cm*. 

Tfie  iron  membranes  (0.025  cm  lliick)  were  obtained  from 
Goodfellow  with  the  following  composition:  Mn,  0.3%  Si,  0.1%;  C, 
<0.08%;  S,  <0.05%;  and  the  balance  is  iron.  They  were  polished 
successively  down  to  600  emery  paper,  degreased  ultrasonically  in 
acetone,  annealed  at  900X  for  2  h  in  pure  hydrogen,  and  furnace 
cooled  In  the  same  atmosphere. 

The  chargirig  solution  was  0.1  A/H^SO,  4-  0.9  A/NajSO^,  pH  3  I 
1.7.  The  solutions  were  prepared  from  analytical  grade  chemicals 
and  double-distilled  water.  They  were  pre-electrolyzed  at  3  mAfor 
24  h  in  external  cells.  All  solutions  were  continuously  deaerated  . 
witfi  pure  hydrogen.  | 


cel!  I 


iron  membrane 

-t  cell  2 


Fig.  1.  Schematic  Illustration  showing  (a,  top)  the  position 
of  the  Iron  membrane  between  the  two  electrolytic  cells  and 
(b,  bottom)  the  boundary  conditions  Imposed  on  the  mem¬ 
brane  by  the  two  cells. 
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The  exit  side  of  the  membrane  was  coated  with  a  thin  layer  of  Pd 
using  a  commercial  electroless  plating  bath.  Pallamerse.  This  Pd 
layer  prolecis  the  exit  side  of  the  membrane  and  provides  a  cat- 
^ylic  surface  for  the  oxidalion  of  hydrogen.  The  exit  side  of  the 
membrane  was  polarized  anodically  to  i  150  mV  (Hq/HoO)  Tlie 

ivTfin  '•  decayed  to  a  constant  value,  normal- 

y  below  0-2  ,iA  cm  .  Tfie  solution  was  introduced  into  tlie  charg- 
ng  s  de  of  the  cell  under  conditions  of  cathodic  protection  to  pro¬ 
tect  the  sample  from  corrosion  if  left  at  the  open-circuit  potential 

JarnrH  T  T  '''®  ex*'  surface)  was 

obtein^H  ®  ®'ee''y-s'a'e  permeation  current  was 

obtained  The  cathodic  current  at  the  charging  surface  was  then 

h  ®  f  reaching  a  steady  permeation 

charninn®p|prf  ‘i®t'  experiments  were  conducted  with 

charging  electrolytes  containing  different  BTA  concentrations. 

Results  and  Discussion 

^  Pe^e^eation  transients  obtained  at  a 

SnPnnTuri'^fRTr'  '  'or  Charging  electrolytes 

®  ^  "e.'^e"'  BTA  concentrations.  The  permeation  current 
ncreases  wi  h  tirne  toward  a  limiting  steady-state  value,  /  The  liq- 

^s®with®!nrIpfcL  f '’ydroger.  through  irondecreal 
k  concentration.  Further,  BTA  shifts  the 
S  "^e  charging  side  of  the  membrane  toward  more 
alhodic  values  to  an  extent  which  also  depends  on  the  BTAcon- 

®®''®''e'  permeation  transients  were 
r  different  charging  currents  and  were  repeated  at  dif- 

&renl  BTA  concentrations. 

The  rate  of  proton  discharge  (reaction  1)  4  is  given  by 


(I  -  «i.(  -^5!) 


|5J 

V  “  *  / 

here  A,  is  the  rate  constant  of  reaction  1,  C„.  is  the  hydrogen  ion 
mcentralion.  ii,,  is  the  degree  of  surface  coverage  of  the  mem- 
arte  with  adsorbed  hydrogen,  »  is  the  transfer  coefficient  of  the 
sct^rge  reaction  (Eq.  1).  -ij  is  the  cathodic  polarization,  Fis  the 
“cstant  flis  the  universal  gas  conslani,  and  Tis  llie 
•solute  temperature.  The  rata  of  llie  (Tale!)  recombination  renc 
■n  (Eq.  2)  IS  given  by 


-  FA.ii?, 


|6) 


For  a  mechanism  of  tfie  HER  involving  coupled  proton  discharge 

^nicuaT®  ®  ''^dous  theoretical  Irenl- 

•nls  relate  the  steady-state  permeation  (;.)  and  recombina- 
1  (/,)  currents  by 


[JT 


[7a] 


ere  D  is  the  diffusivity  of  hydrogen  in  the  membrane  and  L  is  its 
ikness.  In  view  of  Eq.  6,  the  above  equation  also  can  be 
iressed  as 


Equation  7a  predicts  a  straight  line  relation  between  and 
. wtiere  ',  =  4  “  Figure  3  shows  several  such  plots  at  different 
BTA  concentrations.  The  plots  show  satisfactory  straight  lines 
passing  througfi  lire  origin  which  confirm  that  the  reaction  pro- 
ceeds  via  me  Tafel  recombination  and  tliat  the  permeation  rate  is 
controlled  by  hydrogen  diffusion  through  ttie  membrane  and  not  by 
slow  kinetics  at  the  surface,  i.e..  Eq.  4  is  in  local  equilibrium 
Flintier,  tfie  slopes  of  the  straight  lines  decrease,  indicating  lower 
permeation  currents  as  the  BTA  concentration  increases.  Ac- 
fpn?m^  °K  equations,  a  decrease  in  the  permeation  cur¬ 

rent  may  be  brouglit  about  by  a  decrease  in  0„  and/or  in  the  ratio 
".te/k*,.  w'l'?'}  's  the  equilibrium  constant  of  reaction  4.  This  indi- 

w«h  ad^rhpIiVH  '^®®''e®ees ''’e  coverage  ol  the  iron  surface 
wilfi  adsorbed  hydrogen  and/or  diminishes  the  rate  constant  of  the 

ohrp^prt -ffiTA  principle,  a  similar  effect  may  be 

Jrea^fior!  4^  increase  /q„..  i.e.,  to  favor  the  reverse  step 

Figure  4  shows  the  E-log(/)  relations  of  the  HER  obtained  on  iron 
m  the  presence  of  various  concentrations  of  BTA.  Clearlv  BTA 

tXPcP®®  ®  pe'enlial.  It  also  increases  the 

afel  slope  (2.3  FtTUxF)  from  120  tnV  per  current  decade  for  the 
blank  electrolyte  m  146,  169,  and  169  mV  for  the  electrolytes  con- 

hai'"uo^’  mMBTA,  respectively.  These  plots^indicafe 

lhal  tfie  raie-delermining  step  of  ttie  HER  is  a  charge-transfer 
reaction  (r.e.  proton  discharge,  reaction  1)  with  a  transfer  ooelh- 
cient «  equal  to  0.5  lor  tfie  blank  electrolyte  and  «  equal  to  0.35  n 
the  presence  of  10  mMBTA.  Reaction  2  is  excluded  as  a  possible 
rale  determining  step,  as  this  requires  a  Tafel  slope  of  only  30  mV 


FT) 

k.. 


. /i  f *  sieaay-state  hydrogen  permeation  cur- 

(7b]  rent  (I,)  vs.  square  root  of  the  recombination  current  d)  at  dif¬ 
ferent  BTA  concentrations.  ' '' 
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I  plots  of  the  hydrogen  evolution  reaction  at  dif- 

ferent  BTA  concentrations. 


per  current  decade,'  which  was  not  observed  here.  This  analysi' 
o(  Fig,  4  Indicates  that  BTA  also  Inhibits  the  HEP  by  decreaslnq  litr 
rate  constant  ol  proton  discharge  (k,.  see  Eq.  5),  or  the  fraction  o 
uncovered  electrode  area,  which  becomes  (1  -  n„  -  ()„,.)  in  prin. 

r'ilrm, ?  fi' observed  on  a  dect ease  in  the  rale  ol 
recombination  (Eq.  2),  If  llie  mectianisrn  of  the  HER  is  coupled  pro¬ 
ton  discharge-Tafel  recombination  (Eq.  1  and  2)  with  both  steps 
proceeding  at  comparable  rates.  The  Increase  In  Ihe  Talel  slope  is 
determined  solely  by  a  decrease  In  n  (Eq.  5).  ^ 

^  *’’®  between  the  sleady-slnte  perme- 

nine  h  fod  Ihe  caihodic  potential  at  different  BTAconcentra- 
'bo  permeation  current  increases  at  more  neg- 
ninn  ^''e'be'''e)  Polemials  and  decreases  with  Increasing  BTA 
concentration.  This  indicates  that  BTA  Inhibits  Ihe  hydrogen 

the  HER  r/oT  f?  of  ifihiblng 

Jbis  effect  of  BTA  Is  different  than  that  of 
urea  wfiich  inhibits  Ihe  HER  and  promotes  hydrogen 

airi°in inf  promotes  hydrogen  abLrpfon 

®  'b®  blER.  The  mechanism  by  which  such 
rpnoir®  P'?'^®'®  hydrogen  absorption  was  discussed 

immn/  "  o' e'sclronic  interactions  between  Ihe  adsorbed 
p  omolers  and  the  adsorbed  hydrogen  atoms  on  the  metal  surface. 
More  work  must  be  done  to  unravel  tire  mechanism  ol  this  dual 

Ihe  simultaneous  hydrogen  evolution 
and  hydrogen  absorptron  reactions. 


Fig.  5.  Steady-slate  permeation  current  (/) 
potential  (~f)  at  different  BTA  concentrations. 


vs.  cathodic 


Table  (.Values  of  the  steady-state  permeation  current,  ar 
the  concenlrallon  of  hydrogen  In  Iron  at  the  Input  surface 
C\  at  different  concentrations  of  the  BTA  at  a 
potential  of  -0.65  V  (SCE). 


(BTA).  mA4 

L  (it A  cm  ^) 

C  (mol  cm' 

0 

2,17 

4.7  X  10- 

1 

0.46 

1.0  X  lO"' 

10 

0.24 

5.0  X  10  ' 

50 

0.16 

3.5  X  10  ' 

The  steady-stale  permeation  current  is  related  to  the  concentr 
lion  C''  by 

■  t  U/k-r  /  L. 


where  0  IS  Ihe  diffusiviiy  of  hydrogen  in  the  membrane  of  thickne; 
L.  Equation  8  reveals  that  llie  decrease  in  the  steady-slate  perm 
aliori  current  is  brougftl  about  by  a  decrease  in  C".  This  can  I 
caused  by  eillier  a  decrease  in  0„  or  by  a  shift  of  the  equilibria 
portion  of  reaction  4  toward  the  left,  i.e.,  away  from  fowai 

Tins  conclu.sion  is  consistent  with  that  obtained  from  am 
typing  Ihe  data  in  Fig.  3,  in  Ihe  liglil  of  Eq.  7a  and  b. 

Equation  8  can  be  used  to  calculate  Ihe  values  of  C"  under  va 
lotis  conditions  provided  Ihe  value  ol  D  is  known.  We  measured 
dillasivilyof  1.2X  to  *cm’s  '  lor  hydrogen  in  Iron  using  Ihe  equt 
bon  of  breakllirough  lime  given  by  Dovannlhan  and  Slachurski.l 
Table  I  lists  llie  values  of  C"  so  obtained  at  ■  650  mV  (saturate' 
calomel  elecirode  (SCE))  as  a  funclion  of  BTA  concenlrallon. 
the  BTA  concentration  increases,  C”  decreases.  A  10  mM  BtJ 
solulion  decreases  C"  by  about  one  order  of  magnitude.  i 

Conclusion 

The  present  work  sliows  that  BTA  is  an  efficient  inhibitor  for  !h( 
absorption  of  H  into  Iron  and  also  for  tlie  hydrogen  evolution  reac 
ion.  The  iricrease  of  the  BTA  concentration  causes  an  increase  ir 
the  cathodic  polarization  and  In  the  cathodic  Tafel  slope  on  iron.  I 
also  causes  a  decrease  in  Ihe  steady-stale  permeation  current  ant 
hence  n  the  concentration  of  hydrogen  In  the  metallic  lattice.  Thi 
adsorption  of  BTA  on  Hie  Iron  surface  leads  to  a  marked  decrease 
in  IJie  rale  of  proton  discliarge  and  a  shift  of  Hie  surface  equilibril 
urn  M  -  toward  Ihe  left  side  which  results  In  s 

decrease  in  the  concentration  of  absorbed  hydrogen  within  the 
metal  This  effect  is  opposite  to  that  reported  lor  some  othei 
inhibitors  such  as  thiourea  which  Inhibits  Ihe  hydrogen  evolution 
promotes  hydrogen  permeation.  The  specifid 
effect  of  BTA  may  yet  be  sorted  out  from  additional  measuremenid 
and  analysrs.  ! 
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The  Effect  of  Iodide  Ions  on  the  Kinetics  of  Hydrogen 
Absorption  by  Iron 

M.  H.  Abd  Elhamid ,*-»  B.  G.  Ateya ,**’•>  and  H.  W.  Pickering***’^ 

Department  of  Materials  Science  and  Engineering,  The  Pennsylvania  State  University,  University  Park,  Pennsylvania 
16802,  USA 


The  effect  of  iodide  ions  on  the  kinetics  of  both  the  hydrogen  evolution  reaction,  HER,  and  the  hydrogen  absoiption  reaction,  HAR, 
on  iron  (steel)  membranes  was  investigated  using  an  electrochemical  hydrogen  permeation  cell.  Iodide  ions  inhibit  the  HER  by 
decreasing  the  rate  constant  of  proton  discharge.  Iodide  enhances  the  HAR,  and  in  turn  increases  the  concentration  of  hydrogen 
inside  the  metal,  although  it  decreases  the  hydrogen  surface  coverage.  The  increased  HAR  rate  is  due  to  an  increase  in  the  com¬ 
posite  rate  constant  of  the  hydrogen  absoiption  step  at  the  metal  surface. 
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The  presence  of  some  additives  in  the  environment  affects  the 
kinetics  of  the  hydrogen  evolution  reaction  (HER)  on  the  surface 
(see  Eq.  1  and  2)  and,  consequently,  the  hydrogen  absorption  reac¬ 
tion  (HAR)  within  metals  (Eq.  3).  While  some  compounds  inhibit 
both  reactions,  for  example,  benzonitrile,  benzotriazole,  and 
pyrroles, others  inhibit  the  HER  while  promoting  the  HAR, 
e.g.y  iodide  ions,  thiourea,  and  aminotriazole.h2,8,i5,i8-2i  phe  mech¬ 
anism  of  action  of  these  additives  is  attracting  considerable  atten¬ 
tion.  Inhibitors  of  both  the  HER  and  HAR  are  compounds 
which  adsorb  preferentially  on  the  metal  surface  and  hence  decrease 
the  hydrogen  surface  coverage,  which  in  turn  decreases  the  rate  of 
both  the  HER  and  HAR. ^27, 33, 34  gy  similar  argument,  others®'^*^ 
attributed  the  effect  of  promoters  of  the  HAR  (Eq.  3)  to  inhibition  of 
the  recombination  reaction  (Eq.  2).  This  leads  to  a  decrease  in  the 
rate  of  the  HER  (Eq.  2)  but  to  an  increase  in  the  hydrogen  surface 
coverage  and  hence  to  an  increase  in  the  hydrogen  absorption.  It  has 
also  been  suggested  that  the  adsorbed  species,  including  iodide  in 
acid  sulfate  solution,  weaken  the  metal  hydrogen  bondhio»i5,i7 
hence  facilitate  the  entry  of  hydrogen  into  the  metal.  Another  group 
of  authors  attributed  the  accelerating  affect  of  these  additives  on  the 
HAR  to  an  electronic  interaction  between  the  adsorbed  hydrogen 
atoms  and  the  adsorbed  additive While  additives  with  elec¬ 
tron  donating  nature,  such  as  iodide, are  reasoned  to  accelerate 
hydrogen  absoiption,  those  with  electron  withdrawing  nature  have 
the  opposite  effect.^"^’^^'^^ 

Recently,  Qian  et  al?^*^^  investigated  the  effect  of  arsenic  oxide 
and  hydrogen  sulfide  in  promoting  hydrogen  absorption  into  iron 
and  steel  in  alkaline  solutions.  Their  results  revealed  that  arsenic 
oxide  enhances  hydrogen  absorption  through  its  reduction  to  arsenic 
hydride  (ASH3).  These  authors  concluded  that  ASH3  is  desorbed 
from  the  metal  surface,  thereby  increasing  the  number  of  adsorption 
sites  for  hydrogen.  Thus,  the  hydrogen  surface  coverage  and,  conse¬ 
quently,  the  hydrogen  absorption  rate  are  increased^^  as  reported  ear¬ 
lier.*  In  the  case  of  hydrogen  sulfide  in  alkaline  solution,  they  con¬ 
cluded  that  HS”  adsorbs  on  the  surface,  thereby,  decreasing  both  the 
hydrogen  surface  coverage  and  the  hydrogen  desorption  rate,  al¬ 
though  enhancing  hydrogen  absorption.^^  In  a  prior  study  of  the 
Fe/H2S-acid  system,  Iyer  et  alf^  concluded  that  H2S  caused  an 
increase  in  the  hydrogen  surface  coverage  and,  in  turn,  accounted  for 
the  increase  in  the  hydrogen  absoiption  rate. 

The  aim  of  the  present  study  is  to  analyze  the  effects  of  iodide  ions 
on  the  kinetics  of  both  the  HER  and  HAR  on  iron  in  a  more  in-depth 
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manner  than  in  the  past.  Of  particular  importance  is  the  effect  of  iodide 
ions  on  the  hydrogen  coverage  of  the  iron  surface  and  the  rate  constant 
of  the  hydrogen  absorption  step  which  affects  the  rate  of  hydrogen 
permeation  within  iron.  The  analysis  of  the  results  makes  use  of  mod¬ 
els  that  have  been  extensively  tested  in  the  literature. 

Model 

Consider  a  coupled  (Volmer)  discharge-(Tafel)  recombination 
mechanism  for  the  HER,  taking  place  on  the  left  side  of  the  mem¬ 
brane  shown  in  Fig.  1,  i.e. 


H+  +e-  +  M-^M-H,ds 

[1] 

ko 

+  M-H3as  — 2M  +  H2 

[2] 

where  M-H^ds  refers  to  an  adsorbed  hydrogen  atom  on  the  metal  sur¬ 
face  and  k]  and  k2  are  the  rate  constants  of  the  discharge  and  recom- 


Figure  1.  Schematic  illustration  showing  the  boundary  conditions  imposed 
on  the  metal. 
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bination  steps  of  the  HER.  This  adsorbed  hydrogen,  is  also 

involved  in  the  hydrogen  absorption  reaction,  i,e. 

^abs 

Habs+  n  [3] 

*des 


Combining  Eq.  5  and  9  gives  a  relation  between  the  charging 
function  and  the  steady  state  hydrogen  permeation  current 


where  Ha^s  refers  to  the  absorbed  hydrogen  atom  within  the  metal 
lattice  and  Argbs  and  refer  to  the  rate  constants  of  the  forward  and 
backward  directions  of  Reaction  3,  respectively  (see  Fig.  I).  In  addi¬ 
tion  to  the  assumption  of  a  coupled  Volmer-Tafel  mechanism  for  the 
HER,  we  also  assume  the  following:  (i)  that  Reaction  1  (Volmer 
reaction)  proceeds  irreversibly,  (n)  a  Langmuirian  adsoiption  for 
hydrogen,  and  (iii)  hydrogen  diffuses  as  atomic  hydrogen  within  the 
iron  membrane.  These  assumptions  have  been  widely  adopted  in  the 
literature.  As  shown  below,  the  results  of  this  work 

appear  compatible  with  the  above  assumptions.  The  questions  of 
Langmuirian  adsorption  is  addressed  in  the  Conclusion.  According 
to  Eq.  1  the  rate  (current  density)  of  the  HER  is  given  by 

'c  =  -  0h)  [4] 


where  is  the  hydrogen  ion  concentration,  a  is  the  transfer  coef¬ 
ficient  of  the  HER,  and  r\  is  the  hydrogen  overpotential.  According 
to  the  IPZ  analysis,'*^’'^^  Eq.  4  gives 

‘c  -  9h)  [5] 


where  the  left  side  of  Eq.  5  is  referred  to  as  the  charging  function  and 
Zo  =  —  /c/(l  —  0h)  where  i^  is  the  exchange  current  densi¬ 

ty  of  the  HER  and  6^  is  the  hydrogen  surface  coverage  at  equilibri¬ 
um.  The  rate  of  the  recombination  reaction,  Eq.  2,  is  given  by 

»V  =  Fk2el  [6] 


The  recombination  current  density,  4,  is  the  difference  between  the 
cathodic  current  density  and  the  steady  state  hydrogen  permeation 
current  density,  The  steady  state  hydrogen  permeation  current 
density,  z„,  is  given  by 


FDCf 

L 


[7] 


where  D  is  the  diffusion  coefficient  of  atomic  hydrogen  within  the 
metal,  C®  is  the  concentration  of  the  absorbed  hydrogen  just  beneath 
its  surface,  and  L  is  the  thickness  of  the  membrane.  This  current  is 
given  by  the  difference  between  the  rates  of  the  partial  reactions 
involved  in  Eq.  3  1,40,41-43  ^  ^ 


Combining  Eq.  7  and  8  gives 


[8] 


fee  -  [9] 

where  k  is  the  kinetic-diffusion  constant  and  is  defined  as 


k  - 


^abs 


1  + 
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[10] 


This  composite  rate  constant,  k,  in  Eq.  1 0  is  an  improved  version  of 
that  (E')  in  the  original  paper  where  k!'  =  kLID.^^  Equation  10  was 
also  recently  derived  by  Ramasubramanian  et  alf^  Combining  Eq.  6 
and  9  gives 


This  relationship  has  been  experimentally  observed  by  Daftt  et  al} 
and  others  in  various  systems.  Equation  1 1  predicts  a  straight  line  re- 
lation  passing  by  the  origin  between  z»  and  with  a  slope  of  kjF/f^. 


Equation  12  relates  the  charging  function,  z^  exp(F/RT  ari),  which 
includes  the  charge  transfer  kinetics  at  the  charging  surface  of  the 
membrane,  and  the  steady  state  hydrogen  permeation  current  densi¬ 
ty  (Zoo)  which  is  dependant  on  the  transport  properties  of  the  mem¬ 
brane.  It  predicts  a  straight  line  relation  between  the  charging  func¬ 
tion  and  Zoo,  the  slope  of  which  equals  —VJFk  and  the  intercept  is 
Consequently,  values  of  k  and  Zq  can  be  evaluated  graphically  from 
such  plots.  The  value  of  zj,  can  be  used  to  obtain  k^  =  (z!/^^h+) 

^0  while  k  can  be  used  to  evaluate  the  hydrogen  surface  cover¬ 
age,  0^,  from  Eq.  9.  The  value  of  k  can  also  be  used  to  evaluate  the 
recombination  rate  constant,  ^2,  from  the  slope  of  the  results  plotted 
according  to  Eq.  11.  Substituting  for  z^o  from  Eq.  7  and  Eq.  8  and 
rearranging,  one  obtains^  ^43 

C®  =  k^^KD  [13] 

It  follows  that  once  ^  and  0^  are  known,  the  concentration  of  hydro¬ 
gen  inside  the  metal  at  the  charging  side,  C®,  can  be  obtained  from 
Eq.  13,  as  well  as  from  Eq.  7  using  the  measured  and  D  values. 

Experimental 

The  electrochemical  hydrogen  permeation  apparatus  consisted  of 
two  matching  three- electrode  cells,  similar  to  that  used  by 
Frumkin^’  and  by  Devanathan  and  Stachurski.*^®  It  was  used  to  col¬ 
lect  data  on  both  the  HER  and  HAR  on  iron  membranes  in  an  acidic 
solution  of  0.1  N  H2SO4  +  0.9  N  Na2S04  of  pH  1 .8.  The  iron  (steel) 
samples,  0.25  mm  thick,  were  obtained  from  Goodfellow  with  the 
following  composition:  Mn  0.3%;  Si  0.1%;  C  <0.08%;  S  <0.05%; 
and  the  balance  is  iron.  They  were  polished  down  to  0.5  pm  alumi¬ 
na  and  then  degreased  in  acetone  and  washed  with  double-distilled 
water.  They  were  subsequently  annealed  in  pure  hydrogen  at  900®C 
for  2  h  in  a  tube  furnace  and  fhmace  cooled  in  the  same  atmosphere. 
All  solutions  were  prepared  from  analytical  grade  chemicals  and 
double-distilled  water.  Before  admitting  the  solutions  to  the  cell, 
they  were  pre-electrolyzed  at  3  mA  for  2  h  to  remove  impurities  that 
could  otherwise  affect  the  quality  of  the  data.  The  solutions  were 
subsequently  deaerated  with  hydrogen.  The  exit  surface  of  the  mem¬ 
brane  was  coated  with  electrolessly  deposited  palladium  using  Pal- 
lamerse  solution  before  the  membrane  was  mounted  between  the 
two  cells.  The  exit  cell  contained  0.1  N  NaOH  and  the  potential  of 
the  Pd  coated  exit  side  was  set  at  0.150  mV  (Hg/HgO)  using  a  poten- 
tiostat  to  ensure  complete  oxidation  of  the  dissolved  hydrogen  dif¬ 
fusing  through  the  membrane.  In  order  to  minimize  the  production 
of  traps  in  the  membrane  during  the  measurements,  the  cathodic  cur¬ 
rent,  Zp,  at  the  charging  surface  was  set  at  the  highest  value  and  then 
decreased  in  a  setpwise  fashion  after  reaching  within  minutes  the 
steady-state  permeation  current  at  each  value  of  z^.  Similarly,  the 
polarization  curve  measurements  in  the  same  cell,  the  cathodic  cur¬ 
rent,  z’c,  was  decreased  in  a  stepwise  fashion  after  reaching  in  sec¬ 
onds  the  stationary  electrode  potential  of  the  charging  surface  at 
each  value  of  z^.  Further  details  about  the  experimental  setup  are 

reported  elsewhere.^  ^49 

Results  and  Discussion 

Figure  2  shows  the  effect  of  iodide  ion  concentration  on  the  rela¬ 
tion  between  the  steady-state  hydrogen  permeation  current  density, 
z„,  and  the  overpotential,  t],  of  the  HER  at  the  charging  side  of  the 
membrane  (ti  =  E  —  where  E^^  —  “0.349  V  (standard  calomel 
electrode)  for  the  experimental  system  in  this  paper).  Iodide  ions 
enhance  hydrogen  absorption,  z.e.,  the  steady-state  hydrogen  perme¬ 
ation  current  density  increases  for  increasing  iodide  concentration  at 
a  given  oveipotential  of  the  HER.  Alternatively,  at  a  certain  concen¬ 
tration  of  iodide  ions,  larger  negative  overpotentials  produce  larger 
steady-state  permeation  currents,  in  agreement  with  the  results  of 
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Figure  2,  The  relation  between  the  steady-state  hydrogen  permeation  current 
and  the  overpotential  of  the  HER  obtained  on  an  iron  membrane  0.25  mm 
thick  in  0. 1  N  H2SO4  +  0.9  N  Na2S04,  at  different  iodide  ion  concentrations. 


le,  mA  cm'* 

Figure  3.  Tafel  plots  obtained  on  an  iron  membrane  0.25  mm  thick  in  0.1  N 
H2SO4  +  0.9  NNa2S04  at  different  iodide  ion  concentrations. 


Other  authors. The  difhisivity  of  hydrogen  within  the  iron  mem¬ 
brane  was  estimated  using  the  breakthrough  time  method.*^^  A  value 
of  1.2  X  10""^  cm^  s"^  was  obtained  which  is  in  good  agreement  with 
the  values  reported  in  the  literature  for  the  diffusivity  of  hydrogen  in 
iron.^^  The  increase  in  the  steady-state  hydrogen  permeation  current 
density  is  accompanied  by  a  decrease  in  the  rate  of  the  HER,  as 
shown  by  Fig.  3  which  indicates  that  the  iodide  ions  inhibit  the  HER. 

Figure  4  shows  the  relation  between  the  steady-state  hydrogen 
permeation  current  density  and  the  square  root  of  the  hydrogen 
recombination  current  density.  This  figure  shows  data  which  lie  fair¬ 
ly  well  on  the  straight  lines  drawn  through  the  point  of  origin,  in 
accord  with  the  coupled  discharge-recombination  mechanism  of  the 


0  0.2  0.4  0.6  0.8  1  1.2 

Vi,  mA  cm 

Figure  4.  Relationship  between  the  steady-state  hydrogen  permeation  cur¬ 
rent,  /«,  and  the  square  root  of  the  hydrogen  recombination  current, 
obtained  on  an  iron  membrane  in  0.1  N  H2SO4  +  0.9  N  Na2S04,  at  different 
iodide  ion  concentrations. 


HER,^^  in  the  absence  and  in  the  presence  of  iodide  ions.  The  slopes 
of  the  lines  in  Fig.  4,  kjF/J^  (from  Eq.  1 1),  increase  as  the  iodide  ion 
concentration  increases,  see  Table  I.  Thus,  an  increase  in  the  slo^s 
of  these  straight  lines  results  from  an  increase  in  the  ratio  of  k/Jki. 
This  may  be  caused  by  an  increase  in  k  (z.e.,  enhancement  of  the 
HAR,  Eq.  3).  It  may  also  result  from  a  greater  increase  in  k  than  in 
(k2,  and/or  a  decrease  in  ki  (fe.,  retardation  of  the  recombination 
step,  Eq.  2). 

The  Tafel  slope  for  the  blank  solution  (Fig,  3)  was  120  mV  cor¬ 
responding  to  a  transfer  coefficient  a  =  0.5.  This  value  was  used  to 
calculate  the  charging  function,  exp(a'Ti  F/RT)  under  various  con¬ 
ditions.  Figure  5  shows  the  relation  between  the  charging  function 
and  the  steady-state  permeation  current  at  different  iodide  concen¬ 
trations.  The  relations  are  satisfactory  straight  lines  with  negative 
slopes,  —VJFky  in  agreement  with  the  prediction  of  Eq.  12.  The  mag¬ 
nitude  of  the  slope  decreases  shaiply  with  the  iodide  ion  concentra¬ 
tion,  see  Table  I.  The  decrease  in  the  slopes  of  the  lines  can  be 
brought  about  by  an  increase  in  a  decrease  in  Zq,  and/or  a  greater 
increase  in  ^  than  in  Vq.  The  intercepts  of  the  lines  in  Fig.  5  with  the 
y-axis  give  the  values  of  the  exchange  current  density  of  the  HER, 
^0  =  (1  “  ®h)  h  ^^ses  of  05  «  1.  The  discharge  rate  con¬ 
stant  was  then  calculated  from  the  relation  Iq  =  The  slopes 

of  the  lines  in  Fig.  5  were  solved  with  the  intercepts  for  the  values  of 
the  kinetic-diffusion  constant,  ^(according  to  Eq.  12).  The  values  of 
k  at  different  iodide  ion  concentrations  were  then  used  to  estimate 
the  recombination  rate  constant  of  the  HER,  ^2,  from  the  slopes  of 
the  lines  in  Fig.  4,  in  accord  with  Eq.  1 1.  Table  I  shows  values  for  Zq, 
ku  k2.  and  k  at  different  iodide  ion  concentrations.  The  results  show 
that  iodide  ions  inhibit  the  HER  by  decreasing  the  discharge  rate 
constant,  ky  Similar  results  were  found  in  the  literature  for  the  effect 
of  organic  inhibitors  on  the  discharge  rate  constant  of  the  HER.^^»^^ 
On  the  other  hand,  iodide  promotes  the  HAR  by  increasing  the  kinet¬ 
ic-diffusion  constant,  k  (Eq.  9)  with  in  turn  increases  the  concentra¬ 
tion  of  hydrogen,  C®,  just  inside  the  metal  (Eq.  13  and  Fig.  6)  and  the 
rate  of  hydrogen  permeation  (Eq.  7) 

The  results  in  Table  I  also  show  that  the  iodide  ions  slightly 
increase  the  recombination  rate  constant  of  the  HER,  ^2*  The  values 
of  k2  in  Table  I  differ  by  much  more  than  the  standard  deviation  of 


Table  I.  Values  of  4,  kj,  VJFk,  kjf/kj,  and  k  obtained  on  a  iron  membrane  in  0.1  N  H2SO4  +  0.9  N  Na2S04  at  different  iodide  ion  concen¬ 

trations. 


[r],mM 

Zq,  A  cm’’^ 

^j,  cm 

rjFk 

k  [T  cm"' 

pi 

k,  mol  cm"^ 

^2,  rnol  cm“^  s“^ 

0 

2.0  X  10'* 

1.3  X  10"* 

0.40 

1.0  X  10"“' 

5.2  X  10"" 

2.6  X  10"* 

1 

1.2  X  10"* 

7.3  X  10"’ 

0.13 

1.8  X  10"“' 

9.3  X  10"" 

2.7  X  10"* 

10 

9.0  X  10"’ 

5.6  X  10"’ 

0.08 

2.1  X  10"'' 

1.1  X  10"'® 

2.8  X  10"* 

50 

6.0  X  10"’ 

3.8  X  10"’ 

0.03 

3.1  X  10"* 

2.0  X  10"'® 

3.8  X  10"* 

2261 


Journal  of  The  Electrochemical  Society,  147  (6)  2258-2262  (2000) 

80013-4651(99)07-049-4  CCC:  $7.00  ®  The  Electrochemical  Society,  Inc. 


Figure  5.  Relationship  between  the  charging  function,  i^  exp(aFT)//?7),  and 
the  steady-state  hydrogen  permeation  current,  for  an  iron  membrane 
0.25  mm  thick  in  0,1  N  H2SO4  +  0.9  N  Na2S04,  at  different  iodide  ion  con¬ 
centrations. 


n.v 

Figure  6.  Variation  of  the  concentration  of  hydrogen  in  iron  at  the  input  sur¬ 
face,  C°,  with  the  overpotential  of  the  HER  obtained  on  an  iron  membrane 
0.25  mm  thick,  in  0. 1  N  H2SO4  +  0.9  N  Na2S04,  at  different  iodide  ion  con¬ 
centrations. 


the  measured  values  in  three  or  four  duplicate  runs  at  all  P  con¬ 
centrations.  The  mean  values  of  /oo  were  used  to  determine  the  k,  /tj, 
and  k2  values  in  Table  1.  These  changes  in  and  ^2  could  lead  to  a 
decrease  in  the  hydrogen  surface  coverage.  The  values  of  k  were 
used  to  estimate  the  hydrogen  surface  coverage  using  Eq.  9  at  dif¬ 
ferent  values  of  the  hydrogen  overvoltage  and  iodide  concentrations. 
Figure  7  gives  the  relations  between  the  hydrogen  surface  coverage 
and  hydrogen  overpotential.  It  reveals  that  the  hydrogen  coverage 
increases  with  the  increase  in  the  overpotential,  which  confirms 
results  obtained  by  other  authors.^®'^^’^®’"^^"^^’^^’^^  Figure  7  also 
reveals  that  the  hydrogen  surface  coverage  decreases  with  increasing 
concentration  of  iodide  ions  in  the  charging  solution,  in  contrast  to 
the  increase  in  concentration  of  hydrogen  in  the  metal  (Fig.  6). 

Conclusions 

In  acidic  solutions,  iodide  ions  inhibit  the  HER  on  iron  while 
enhancing  the  HAR.  Using  the  IPZ  analysis  of  the  experimental 
results,  we  show  that  iodide  ions  lead  to  the  following  effects. 

1.  A  decrease  in  the  hydrogen  surface  coverage,  0^; 

2.  An  increase  in  the  concentration  of  hydrogen  in  the  iron,  C°; 

3.  A  decrease  in  the  discharge  rate  constant  of  the  HER,  k^,  and 
hence  in  the  exchange  current  density, 

4.  An  increase  in  the  value  of  the  kinetic-diffusion  constant  of  the 
HAR,  ky  which  causes  an  increase  in  the  rate  of  hydrogen  absorption 
within  the  metal  (Eq.  3)  and  in  C®  (Eq.  13). 


Figure  7.  Variation  of  the  hydrogen  surface  coverage,  0h,  with  the  overpo¬ 
tential  of  the  HER  obtained  on  an  iron  membrane  0.25  mm  thick  in  O.I  N 
H2SO4  +  0.9  N  Na2S04,  at  different  iodide  ion  concentrations. 


5.  A  slight  increase  in  the  recombination  rate  constant  of  the 
HER, 

Accordingly,  the  above  data  and  IPZ  analysis  reveal  that  the 
increase  in  the  hydrogen  absorption  rate  in  the  presence  of  iodide 
ions  is  caused  by  an  increase  in  the  kinetic-di^sion  constant,  L 
This  increase  in  k  dominates  over  the  decrease  in  surface  coverage 
that  is  also  caused  by  the  iodide  ion.  Note  that  is  a  composite  quan¬ 
tity  (Eq.  10)  which  includes  the  rate  constants  of  both  the  absorption 
and  desorption  steps  (forward  and  reverse  directions  of  Reaction  3), 
the  difflisivity  of  hydrogen  within  the  membrane  and  the  membrane 
thickness.  Thus,  an  increase  in  k  is  caused  by  increases  in 
jtjjgs  within  the  thermodynamic  constraint  that  their  ratio  remains 
unchanged.  The  extent  of  the  contributions  of  the  changes  in 
fcjjgs  to  the  measured  increase  in  k,  can  be  assessed  based  upon  the 
measurements  of  /«>  for  membranes  of  different  thickness.^3, 44.49, 53 
These  experiments  have  yet  to  be  done  for  iodide  ions. 

This  effect  of  iodide  ions  may  result  from  its  interaction  with  the 
metal  surface  and/or  with  the  adsorbed  hydrogen  atoms.  Further 
measurements  are  needed  to  reveal  the  nature  and  extent  of  this  inter¬ 
action  with  the  metal  surface.  The  interaction  between  iodide  ions 
and  adsorbed  hydrogen  atoms  could  also  be  treated  within  the  con¬ 
text  of  Frumkin  adsoiption.  However,  the  above  changes  in  the  kinet¬ 
ic  parameters  are  only  small  (A2)  or  moderate  (Zq,  /tj,  and  k)  and  are 
not  associated  with  any  significant  deviation  from  the  behavior  of  the 
system  in  the  absence  of  iodide  ions.  Furthermore,  the  different  sets 
of  data  are  internally  consistent  and  compatible  with  the  underlying 
model,  7,30,40-46  which  has  been  widely  adopted  in  the  literature, 
in  the  presence  or  absence  of  additives.  Consequently,  one  might 
argue  in  favor  of  small,  if  any,  lateral  interactive  frumkin)  effects. 
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List  of  Symbols 

Ch-h  hydrogen  ion  concentration,  mol  cm“^ 

C°  concentration  of  the  absorbed  hydrogen  at  the  input  surface,  mol  cm  ^ 
D  hydrogen  diffusion  coefficient,  cm^  s*"’ 

F  Faraday’s  constant,  C  moP^ 

I'i  =  /-i,  Ch+  =  ij{\  -  eg),  A  cm-^ 

/q  exchange  current  density,  A  cm*”^ 

Zj.  steady-state  recombination  current  density,  A  cm"^ 

steady-state  hydrogen  permeation  current  density,  A  cm“^ 

/tj  discharge  rate  constant,  cm  s“  ^ 
k2  recombination  rate  constant,  mol  cm“^  s"  * 
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/TjjIjj  absorption  rate  constant,  mol  cm“^ 

^dcs  desorption  rate  constant,  cm 

k  kinetic-diffiision  constant,  mol  cm~^  s”  • 

R  general  gas  constant,  8.3 14  J  moP  ‘  K”  • 

T  absolute  temperature,  K 

Greek 

a  transfer  coefficient,  dimensionless 

7)  cathodic  overvoltage,  V 

hydrogen  surface  coverage  at  equilibrium,  dimensionless 
0H  hydrogen  surface  coverage,  dimensionless 
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Calculation  of  the  Hydrogen  Surface  Coverage  and  Rate  Constants  of  the 
Hydrogen  Evolution  Reaction  from  Polarization  Data 
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We  show  in  this  communication  that  the  polarization  data  of  the  hydrogen  evolution  reaction  (HER)  can  be  analyzed  to  calculate 
the  hydrogen  surface  coverage  and  the  rate  constants  of  the  hydrogen  discharge  and  recombination  reactions  for  metals  which  have 
very  low  permeabilities  of  hydrogen  and  on  which  the  HER  proceeds  through  a  coupled  Volmer  discharge-Tafel  recombination 
mechanism.  The  analysis  is  applied  to  the  results  of  the  HER  on  copper  and  iron.  For  both  metals,  this  polarization  analysis  yields 
exchange  current  densities  that  are  comparable  to  the  literature  values  and  degrees  of  surface  coverage  which  vary  with  potential 
in  a  manner  that  is  well  documented  in  the  literature.  The  results  on  iron  yield  rate  constants  that  are  in  a  good  agreement  with 
those  predicted  in  other  studies  using  the  Iyer,  Pickering,  and  Zamanzadeh  analysis. 
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The  hydrogen  evolution  reaction  (HER)  is  of  particular  impor¬ 
tance  in  electrochemistry  in  view  of  its  relevance  to  corrosion,  elec¬ 
trolytic  hydrogen  production,  electroplating,  electrocatalysis,  etc.^“^^ 
Adsorbed  hydrogen  (Hg^g)  is  an  important  intermediate  in  the  evolu¬ 
tion  of  electrolytic  hydrogen  on  metal  surfaces,^’^*^*"^*^’^^’’'’'^^  and 
also  because  of  its  absorption  within  the  metals  and  the  subsequent 
problem  of  hydrogen  embrittlement.^‘^»^’^^‘*^  Proper  characteriza¬ 
tion  of  the  mechanism  of  the  HER  requires  knowledge  of  the  values 
of  the  rate  constants  of  the  various  steps  and  the  degree  of  surface 
coverage  of  the  metal  with  adsorbed  hydrogen.  The  experimental 
determination  of  the  degree  of  coverage  of  metals  with  adsorbed  hy¬ 
drogen  is  a  task  that  is  fraught  with  many  difficulties.^* 9 

In  this  communication,  we  show  that  the  polarization  curves  of 
the  HER  can  be  analyzed  to  yield  the  individual  rate  constants  of  the 
two  steps  of  the  HER  and  the  degree  of  coverage  of  the  metal  sur¬ 
face  with  adsorbed  hydrogen.  The  analysis  is  suitable  for  metals 
which  have  low  absorption  rates  for  hydrogen  and  for  which  the 
HER  occurs  by  the  coupled  volmer  discharge-Tafel  recombination 
mechanism.  The  proposed  analysis  is  applied  to  the  polarization 
curve  of  the  HER  on  Cu  and  Fe.  Using  iron  under  conditions  of  a  rel¬ 
atively  small  rate  of  hydrogen  absorption,  a  comparison  is  made 
between  results  obtained  using  the  present  polarization  analysis  and 
those  obtained  from  steady-state  hydrogen  permeation  data  using  the 
Iyer,  Pickering,  and  Zamanzadeh  (IPZ)  analysis.^^*^’ 

Experimental 

Cathodic  polarization  experiments  were  carried  out  on  copper 
samples  of  purity  99.5%  and  iron  samples  of  purity  99.5%  in  0.1  N 
H2SO4  +  0.9  N  Na2S04  (pH  1.8),  using  an  EG&G  PAR  potentiostat 
model  276.  The  samples  were  polished  down  to  0.5  \xm  alumina, 
rinsed  with  acetone  and  double  distilled  water.  The  test  solutions 
were  prepared  from  analytical  grade  chemicals  and  double  distilled 
water.  All  solutions  were  pre-electrolyzed  at  3  mA  for  2  h  before 
measurements  were  taken.  Before  running  the  experiments  the  solu¬ 
tions  were  deaerated  using  hydrogen  to  remove  the  dissolved  oxygen 
from  the  solution.  A  conventional  three-electrode  cell  was  used  with 
saturated  calomel  as  a  reference  electrode  and  two  identical  graphite 
rods  as  the  counter  electrodes. 

Results  and  Discussion 

Theoretical  background. — ^Assume  the  HER  proceeds  by  a  dis¬ 
charge-recombination  mechanism,  Le, 
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+  e“  +  M  M-Hgds 

[1] 

kn 

U  +  M-H,d,->2M  +  H2 

[2] 

where  refers  to  an  adsorbed  hydrogen  atom  on  the  metal  sur¬ 

face  and  ki  and  ^2  21*®  rate  constants  of  the  discharge  and  recom¬ 
bination  steps  of  the  HER,  respectively.  This  adsorbed  hydrogen  is 
also  involved  in  an  absorption  reaction,  i.e. 

M-Hgds=^  .  [3] 

where  Hg^s  refers  to  the  hydrogen  absorbed  within  the  lattice  of  the 
metal.  A  measure  of  a  metal’s  tendency  for  absorbing  hydrogen  is  its 
steady-state  rate  of  hydrogen  permeation,  too.  However,  the  i^^  value 
could  underestimate  the  absorption  rate  that  occurs  prior  to  steady- 
state  pemieation,  since  the  transient  rate  of  hydrogen  absoiption  can  be 
much  higher  than  the  steady-state  value,  especially  initially  when  there 
can  exist  a  high  density  of  unfilled  traps  for  hydrogen  in  the  metal,  eg., 
in  the  case  of  certain  steels  or  in  the  case  of  metals  which  exhibit  high 
solubilities  and/or  diffusivities  for  hydrogen,  eg.,  palladium.  Thus,  if 
4  data  are  used  as  a  measure  of  hydrogen  absoiption,  it  is  advisable 
that  4  is  orders  of  magnitude  below  4  in  order  that  the  absorption  rate 
is  insignificant  at  all  times  with  respect  to  the  rate  of  the  HER. 

At  steady  state  the  rate  of  generation  of  by  Reaction  1  (4) 
equals  the  sum  of  the  rates  of  Reactions  2  (i^)  and  3  (4), 


in  ir  "h  ^00 


[4] 


For  some  metals,  such  as  the  face-centered  cubic  metals  copper  and 
nickel,  the  permeabilities  of  hydrogen  are  typically  very  low  in  view 
of  their  low  diffusivities  for  hydrogen  in  comparison  to  the  rate  of 
the  HER,  i.e.,  4  =  4  »  4.  This  condition  (4  s  4)  could  also  hold 
for  other  metals  like  iron  where  the  value  of  4,  though  readily  mea¬ 
surable,  nevertheless  is  often  some  two  to  three  orders  of  magnitude 
smaller  than  4.  However,  for  iron  in  the  presence  of  poisons  such  as 
hydrogen  sulfide,  the  measured  values  of  4  are  much  higher  and  the 
assumption  4  =  4  becomes  less  valid.  Accordingly,  if  4  is  insignif¬ 
icant  with  respect  to  4j  4bs  insignificant  with  respect  to 

4  during  measurement  of  the  steady-state  polarization  curve,  in 
which  case  the  analysis  can  yield  the  hydrogen  surface  coverage  and 
the  rate  constants  of  the  HER,  as  is  shown  below. 

The  rate  of  proton  discharge  (denoted  4)  in  the  Tafel  region  of  the 
polarization  curve  is  given  by 


»c  =  -  9„)  exp|^-^a-nj 

=  /o(l  -  e„)  exp^-- ^a-nj 


[5] 


Journal  of  The  Electrochemical  Society,  147  (6)  2148-2150  (2000) 
$0013-4651(99)10-058*2  CCC:  S7.00  ©  The  Electrochemical  Society,  Inc. 


2149 


Figure  1.  The  relation  between  the  cathodic  current  and  the  electrode  poten¬ 
tial  obtained  on  copper  in  0.1  N  H2SO4  +  0.9  N  Na2S04. 


where  -  9^),  is  the  exchange  current  densi¬ 

ty  of  the  HER  at  the  limit  of  the  equilibrium  coverage,  9^  =  0,  Ch+ 
is  the  hydrogen  ion  concentration,  a  is  the  transfer  coefficient  of  the 
HER,  and  is  the  overpotential  driving  the  HER.  The  rate  of  the 
recombination  reaction,  Eq.  2,  is  given  by 

h  ~  [6] 

Rearranging  Eq.  5  gives^^*^^ 


'c  =  ‘oi^  -  ®H 

:harg 
ER. ! 

ft  =  JiL 

Inserting  Eq.  4  with  «  4,  Eq.  8a  becomes 

o  =  VL 


[7] 


The  left  side  of  Eq.  7  is  called  the  charging  function  which  combines 
the  rate  and  driving  force  of  the  HER.  Solving  Eq.  6  for  9^  gives 


[8a] 


Combining  Eq.  7  and  8b  yields 


ic  exp 


[8b] 


[9] 


Equation  9  predicts  a  straight  line  relation  between  the  charging 
function,  4  Qxp(F/RT arj),  and  From  the  slope  and  the 

intercept  (Q  of  this  plot  one  obtains  i^  and  ^2- 
From  Eq.  5  at  the  limit  of  9h  =  9fi  one  has 

-  F/:,Ch+(1  -  9-h) 


=  F’A:,Ch+,  for  9^  =  0  ,  [10] 

While  4  can  be  inserted  into  Eq.  10  to  calculate  kj,  ^2  can  be  used  to 
calculate  the  hydrogen  surface  coverage  at  different  potentials,  using 
Eq.  6,  assuming  4  =  i^. 


Data  analysis. — ^Figure  1  shows  the  Tafel  plot  obtained  for  hydro¬ 
gen  evolution  on  copper  in  0. 1  N  H2SO4  +  0.9  N  Na2S04  at  25°C.  The 
Tafel  line  has  a  slope  of  121  mV  {i.e.,  a.  ^  0.5)  which  is  comparable 


0.40  0.60  0.80  1.00  1.20  1.40  1.60 


(nnA 

Figure  2.  The  relation  between  the  charging  function  and  the  square  root  of 
the  cathodic  current  obtained  on  copper  in  0.1  N  H2SO4  +  0.9  N  Na2S04. 

to  the  values  in  the  literature.^^’^^  Figure  2  is  a  plot  of  the  charging 
function,  4  txp(F/RT  a-q),  v.y.  It  shows  a  satisfactory  straight  line 
with  a  negative  slope  and  a  positive  intercept  (see  Eq.  9).  Table  I  shows 
the  values  of  the  exchange  current  density,  4,  discharge  rate  constant, 
/T],  and  the  recombination  rate  constant,  k2,  which  were  calculated 
using  the  slope  and  the  intercept  of  the  line  in  Fig.  2.  The  values  are 
4  =  8  X  10“®Acm“^^j  =  5  X  10”®  cms"’^  andA2  =  1.8  X  10““^ 
mol  cm"2  s”^  The  value  of  the  recombination  rate  constant,  was 
used  to  calculate  the  hydrogen  surface  coverage  9^  using  Eq.  6.  Fig¬ 
ure  3  shows  the  relation  between  9^  and  the  electrode  potential.  This 
figure  shows  that  the  hydrogen  surface  coverage  increases  as  the 
potential  becomes  less  noble,  z.e.,  as  the  overpotential  for  the  HER 
increases.  This  trend  is  well  documented  in  the  literature.^^*^^^"* 

The  present  polarization  analysis  was  also  applied  to  the  polariza¬ 
tion  curve  of  hydrogen  evolution  obtained  on  iron  in  deaerated  0.1  N 
H2SO4  +  0.9  N  Na2S04.  Figure  4  shows  the  relation  between  the 
charging  function  and  the  square  root  of  the  cathodic  charging  current 
obtained  on  iron.  The  charging  function  was  calculated  using  a  value 
of  a  ~  0.5  as  measured  in  this  work.  This  figure  shows  a  straight  line 
relation  with  a  negative  slope  in  accord  with  Eq.  9.  A  comparison  is 
made  in  Table  I  between  the  values  of  the  exchange  current  density 
and  the  rate  constants  obtained  from  this  analysis  and  those  obtained 
for  this  same  system  from  steady-state  hydrogen  permeation  data 
using  the  IPZ  analysis.^^*^^  The  table  shows  that  the  present  polariza¬ 
tion  analysis  gives  the  same  order-of-magnitude  values  as  the  IPZ 
analysis  gives.  The  values  of  the  hydrogen  surface  coverage  estimat¬ 
ed  from  the  IPZ  analysis  were  also  compared  with  those  obtained 
using  the  present  analysis,  see  Fig.  5.  This  figure  shows  a  good  agree¬ 
ment  in  particular  at  the  less  noble  potentials.  At  the  more  noble 
potentials,  eg.,  -0.63  V,  the  value  of  9^  calculated  in  this  work 
is  approximately  15%  lower  than  that  obtained  from  the  IPZ  analysis. 

The  results  presented  herein  show  how  the  polarization  data  can 
be  analyzed  for  the  kinetic  parameters  (rate  constants)  and  the  hy¬ 
drogen  surface  coverage  of  the  HER  on  metals,  when  the  hydrogen 
absoiption  rate  is  negligible  with  respect  to  the  cathodic  hydrogen 
evolution  current  density,  4.  Analysis  of  the  results  on  copper  gives 
a  value  for  the  exchange  current  density  of  the  HER  which  is  simi¬ 
lar  to  those  reported  in  the  literature.22.37, 38  jt  was  shown  that  the  hy¬ 
drogen  surface  coverage  changes  with  the  potential  in  a  trend  that  is 
well  documented  in  the  literature.  Applying  the  above  analysis  to  the 
polarization  data  for  iron  gives  values  of  the  exchange  current  den¬ 
sity,  rate  constants,  and  values  for  the  hydrogen  surface  coverage 
which  are  comparable  to  those  obtained  by  applying  the  IPZ  analy¬ 
sis  to  hydrogen  permeation  data. 


Table  I.  Comparison  of  the  exchange  current  density,  discharge  and  recombination  rate  constants,  and  the  transfer  coefficient  for  iron  in  0.1 
N  H2SO4  +  0.9  N  Na2S04  obtained  from  the  IPZ  analysis  and  the  above  polarization  data  analysis.  The  same  quantities  were  also 
obtained  for  copper  using  the  polarization  data  analysis. 

_ 4  (^  cm "2)  /:2  (mol  s“*)  ^i(cms“*)  a 

Analysis  Fe  Cu  Fe  Cu  Fe  Cu  "Pe  ci 


IPZ  2.0  X  10”®  —  2.6  X  10”®  —  1.3  X  10"®  ~ 

Polarization  data  2.0  X  10"®  8.0  X  10"®  3,0  X  10"®  1.8  X  10"’  1.3  X  10“®  5.0  X  10"® 


2150 


Journal  of  The  Electrochemical  Society,  147  (6)  2148-2150  (2000) 
SOOl 3-465 1(99)1 0-058*2  CCC:  $7.00  ©The  Electrochemical  Society,  Inc. 


^.75  -0.8  -0.85  -0.9  *0.95 
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Figure  3.  The  relation  between  the  hydrogen  surface  coverage  and  the  elec¬ 
trode  potential  obtained  on  copper  in  0.1  N  H2SO4  +  0.9  N  Na2S04. 


Figure  4.  The  relation  between  the  charging  function  and  the  square  root  of 
the  cathodic  current  obtained  on  iron  in  0.1  N  H2SO4  +  0.9  N  Na2S04. 


The  above  described  analysis  of  polarization  data  is  clearly  less 
applicable  in  cases  where  hydrogen  absorption  and  permeation  are 
significant  with  respect  to  i^^  e.g.,  in  the  case  of  Pd  which  has  a  very 
high  permeability  for  hydrogen,  and  iron  in  poisoned  solutions  where 
the  presence  of  the  poison  significantly  increases  the  rate  of  hydrogen 
absorption.  The  above  analysis  of  polarization  data  invovles  a  limit¬ 
ing  case  of  the  IPZ  analysis, where  is  very  small  compared  to 
i^.  A  comparison  between  the  results  obtained  on  iron  using  the  IPZ 
analysis  and  the  above  polarization  data  analysis  shows  that  both 
analyses  give  the  same  order-of-magnitude  values  of  the  kinetic  para¬ 
meters  (exchange  current  density  and  recombination  rate  constant) 
and  of  the  hydrogen  surface  coverage.  Being  less  than  the  complete 
IPZ  analysis,  the  above  polarization  data  analysis  cannot  predict  the 
rate  constants  of  the  hydrogen  absorption  and  desorption  reactions. 
Acknowledgment 

The  authors  acknowledge  financial  support  of  this  work  by  the 
U.S.  Steel  Corporation,  the  Office  of  Naval  Research  (grant  no. 
N00014-96-1-0913),  and  the  Division  of  the  International  Program 
of  the  National  Science  Foundation  (grant  INT-9724698). 

The  Pennsylvania  State  University  assisted  in  meeting  the  publication 
costs  of  this  article. 

List  of  Symbols 

Cp{+  hydrogen  ion  concentration,  mol  cm“^ 

F  Faraday’s  constant,  C  moP  * 

=  /■A,C„+  =  iV(l  -  f  „),  A  cm-2 
4  exchange  current  density,  A  cm  ^ 

steady-state  recombination  current  density,  A  cm”^ 

Zp  steady-state  hydrogen  permeation  current  density,  A  cm“^ 
ki  discharge  rate  constant,  cm  s'*  * 
k2  recombination  rate  constant,  mol  cm”^  s~^ 

R  the  general  gas  constant,  8.3 14  J  moP’ 

T  absolute  temperature,  K 

Greek 

ot  transfer  coefficient,  dimensionless 

T)  cathodic  overvoltage,  V 

0^  hydrogen  surface  coverage  at  equilibrium,  dimensionless 

0H  hydrogen  surface  coverage,  dimensionless 
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Determination  of  the  Rate  Constants  of  Hydrogen 
Absorption  into  Metals 
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Hydrogen  absorption  into  iron  (steel)  membranes  was  studied  using  an  electrochemical  hydrogen  permeation  technique.  The  Iyer, 
Pickering,  and  Zamanzaden  analysis  of  the  measured  steady-state  hydrogen  permeation  rates,  obtained  for  different  membrane 
thiclmesses,  gives  the  rate  constants  of  hydrogen  absorption,  Agbs,  and  desorption,  addition  to  those  quantities  which  can  be 

obtained  from  data  for  a  single  membrane  thickness:  z.e.,  the  hydrogen  surface  coverage,  6^,  exchange  current  density,  4,  rate  con¬ 
stants,  ki  and  ^2,  of  the  hydrogen  evolution  reaction,  and  the  important  composite  quantity  called  the  kinetic-diffusion  constant, 
^  =  ^abs^(^  ^des  ^/^)»  where  L  is  the  membrane  thickness  and  D  is  the  diffusivity  of  hydrogen  within  the  metal.  The  product  of 
k  and  0^  determines  the  steady-state  permeability  of  hydrogen  within  the  metal.  For  large  D  and/or  small  L,k-  A-bs,  and  the  hydro¬ 
gen  permeation  process  is  under  surface  control.  From  the  analysis,  we  obtained  k^^^  =  2  X  10“^®  mol  cm“^  s“'  and  -  1.9  X 
10“^  cm  s“^  for  iron  in  acidic  sulfate  solution. 
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Manuscript  submitted  September  20,  1999;  revised  manuscript  received  May  I,  2000. 


The  absoiption  of  electrolytic  hydrogen  into  metals  is  the  first 
and  necessary  step  in  the  hydrogen  embrittlement  of  metals.  It  is 
equally  important  in  the  permeation  and  storage  of  hydrogen  within 
metals.  ^3-19  process  of  hydrogen  absorption  has  long  been  stud¬ 
ied  using  an  electrochemical  permeation  cell.^®'^^  In  this  technique, 
the  hydrogen  evolution  reaction  (HER)  occurs  on  one  side  (charging 
side)  of  a  thin  metal  membrane  (see  Fig.  1) 

.  ki 

H+  +  e’  +  M  M-Hads  [1] 

where  M-H^ds  refers  to  an  adsorbed  hydrogen  atom  on  the  metal  sur¬ 
face  and  ki  is  the  rate  constant  of  the  discharge  step  of  the  HER  This 
adsorbed  hydrogen  atom  combines  with  another  adsorbed  hydrogen 
atom  to  form  a  hydrogen  molecule,  i.e. 


h 

M-H,ds  +  M-H,ds - ^  2M  +  H2(g)  [2] 

where  lo^  is  the  rate  constant  of  the  recombination  step  of  the  HER. 

This  adsorbed  hydrogen  on  the  charging  surface  of  the  membrane 
is  also  involved  in  the  hydrogen  absorption  reaction  (HAR)  by  which 
the  adsorbed  hydrogen  is  absorbed  in  the  metal  lattice  beneath'  the 
surface,  i.e. 


^abs 

M-H3ds-^H3b,  +  M  [3] 

^des  ^  ^ 

where  H^bs  refers  to  the  absorbed  hydrogen  atom  within  the  metal 
lattice  and  and  the  rate  constants^of  the  forward  and 

backward  directions  of  Reaction  3,  respectively  (see  Fig.  1).  This 
absorbed  hydrogen  diffuses  (permeates)  within  the  metal  lattice  and 
is  oxidized  at  the  other  (exit)  side  of  the  membrane,  under  the  effect 
of  an  anodic  potential,  ^ 

The  steady-state  anodic  current  at  the  exit  surface  of  the  mem¬ 
brane  is  a  measure  of  the  steady-state  permeation  rate  of  hydrogen 
through  the  membrane,  i.e. 
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where  D  is  the  diffusion  coefficient  of  atomic  hydrogen  within  the 
metal,  C°  is  the  concentration  of  the  absorbed  hyi-ogen  just  beneath 
the  metal  surface  on  the  charging  side  where  the  HER  takes  place, 
and  L  is  the  thickness  of  the  membrane.  Equation  4  holds  for  hydro¬ 
gen  permeation  when  trapping  effects  are  insignificant,^^  as  in  the 
case  of  the  carefully  annealed  iron  (steel)  used  here.  Much  work  has 
been  published  on  hydrogen  permeation  using  this  technique.^^’^^ 
Analysis  of  the  results  in  most  of  these  works^*^  0,24,26-33  considered 
the  diffusion  of  hydrogen  within  the  metal  membrane  to  be  the  slow 
step  and  the  surface  reaction  (Reaction  3)  to  be  in  quasi  equilibrium. 
Within  this  framework,  the  individual  rate  constants  of  the  forward 
and  backward  directions  of  Reaction  3  could  not  be  evaluated.  Some 
workers"^*^^’^^’^®  were  able  to  take  into  account  the  surface  reaction 
at  the  entry  side  of  the  membrane  (see  Fig.  1)  to  evaluate  the  de¬ 
sorption  rate  constant  (the  reverse  direction  of  Reaction  3).  Howev¬ 
er,  evaluation  of  the  absoiption  rate  constant  (the  forward  direction 
of  Reaction  3)  was  not  possible  without  the  use  of  independent 
methods  to  evaluate  the  hydrogen  surface  coverage.23,30  Within  the 
framework  of  the  Iyer,  Pickering,  and  Zamanzaden  (IPZ)  analy- 
sis,?^’^^  which  does  not  require  that  Reaction  3  be  in  quasi  equilibri- 


C°  -  0  \ 

0H  -  0  /  at  X  =  L 


Figure  1.  A  schematic  illustration  showing  the  conditions  on  the  hydrogen 
entry  (charging)  and  exit  sides  of  the  iron  (steel)  membrane. 
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um,  the  fractional  degree  of  coverage  as  well  as  both  the  absorption 
and  desorption  rate  constants  can  now  be  determined. 

This  paper  presents  the  first  results  of  experimental  measure¬ 
ments  of  the  rate  constants  of  the  absorption  and  desorption  steps  of 
hydrogen  evolution,  using  the  IPZ  analysis,  on  iron  membranes  of 
different  thicknesses  in  a  buffered  acid  medium.  The  analysis  of  the 
experimental  results  also  yields  the  degree  of  hydrogen  surface  cov¬ 
erage,  the  concentration  of  hydrogen  in  the  membrane  at  the  input 
surface,  and  the  extent  of  surface  control  during  the  absorption-per¬ 
meation  process.  We  also  revise  one  of  the  important  equations  of 
the  IPZ  analysis  and  the  conclusions  thereof 

Theory 

According  to  Reaction  1,  the  rate  (current  density)  of  the  HER 
for  Tafel  behavior  is  given  by 


where  is  the  rate  constant  for  proton  discharge,  C„+  is  the  hydro¬ 
gen  ion  concentration,  0^  is  the  fractional  hydrogen  surface  cover¬ 
age,  a  is  the  transfer  coefficient  of  the  HER,  which  can  be  readily 
calculated  from  the  measured  Tafel  slope,  ti  is  the  hydrogen  overpo¬ 
tential,  andi%R,  and  7  have  their  usual  meanings.  Equation  5  can  be 
rearranged  to  give^®*^^ 


where  —  ij{\  —  0g),  is  the  exchange  current  densi¬ 

ty  of  the  HER,  and  0^  is  the  equilibrium  coverage.  The  left  side  of 
Eq.  6  is  referred  to  as  the  charging  function,  which  previously 
appeared  as  a  key  function  in  the  IPZ  analyses  of  the  rate  constants 
of  the  HER  and  HAR  for  iron  and  other  metals  which  readily  absorb 
hydrogen/^-2^ 

The  rate  of  the  recombination  reaction.  Reaction  2,  is  given  by 

h  =  [7] 

The  adsorbed  hydrogen  resulting  from  Reaction  1  is  distributed  be¬ 
tween  Reactions  2  and  3.  Hence,  at  steady  state 

4  ■“  h  ^oo  [8] 

^  charging  side  of  the  membrane  reveals 

that3A21, 28-30 


Combining  and  rearranging  Eq.  4  and  9  gives^*"**^^’^® 


^00  fq 


1  +  h 


'dcs 


D) 


9„ 


[9J 


[10] 


To  sum  up  the  kinetics  and  transport  properties  of  the  system  in  one 
constant,  Eq.  10  is  expressed  as 


ice 

where  kis  a  kinetic-diffusion  constant 

=  ^bs 


[11] 

[12] 


Equation  12  reveals  that  k  (and  hence  see  Eq.  10)  increases  as  L 
decreases  and/or  D  increases.  This  equation  conveys  the  correct 
physical  picture  of  the  dependence  of  k  (and  hence  /„)  on  both  D  and 
Z,  and  is  a  modification  pointed  out  by  Ramasubramanian  et  al?^  of 
the  quantity  (k")  given  in  the  original  IPZ  analysis,  where  k!*  = 
28,29  phe  constant  k  is  important  in  that  it  combines  the  kinet¬ 
ics  of  hydrogen  absoiption  and  desorption  at  the  charging  side  of  the 
membrane  surface  and  the  transport  properties  within  the  membrane 


and  hence  is  called  the  kinetic-diffusion  constant.  Previously,  we 
have  called  k'  the  thickness-dependent  absoiption-desoiption  con¬ 
stant,  but  in  addition  to  its  cumbersome  length,  this  definition 
ignores  the  dependence  of  ^  on  D  which  is  just  as  important  as  that 
on  L  Henceforth,  k  is  referred  to  as  the  kinetic-diffusion  constant. 

At  the  limit  of  large  D  and  small  Z,  such  that  1  »  k  = 

/Tabs,  and  hence  using  Eq.  11,  one  obtains 

“  '^‘^abs^H  [13] 

This  indicates  that  the  rate  of  hydrogen  permeation  is  controlled 
by  the  surface  step.^  However,  at  the  limit  of  large  Z  and/or  small  Z, 
Eq.  12  gives  k  =  Z)tab/^deA  which  when  substituted  in  Eq.  11,' 
yields 


Ico  = 


=  ^abs  Q 


Z  k 


’des 


[14] 


Under  this  condition,  the  absorption-permeation  process  approaches 
pure  diffusion  control,  in  which  case  C®  is  independent  of  Z,  which 
is  a  limiting  case  of  Eq.  4  when  Eq.  3  is  in  quasi  equilibrium. 

The  reciprocal  of  k  (Eq.  12)  is  given  by 


^  ^abs  ^abs  ^ 


[15] 


which  is  a  corrected  version  of  that  given  previously.^®’^^  On  the 
right  side  of  Eq.  15,  the  contribution  of  the  first  term  to  the  sum  of 
the  two  terms  is  a  measure  of  the  amount  of  surface  control  during 
the  absorption-permeation  process.  Once  k  is  known  (see  below), 
one  can  plot  Hk  vs,  the  membrane  thickness,  Z,  to  obtain  a  straight 
line  with  a  slope  of  W-^^abs  and  an  intercept  of  Uk^^^.  The  slope 
and  intercept  can  be  solved  for  k^^^  and  Combining  Eq.  7 

and  1 1  gives^^’^^ 


[16] 

This  relation  has  been  experimentally  observed  by  Dafft  et  al?^  and 
others  in  various  systems.  Equation  16  predicts  a  straight-line  rela- 
tionbetween  and  74  passing  by  the  point  of  origin  with  a  slope  of 
k^FIk^.  The  slope  can  be  used  to  obtain  ^  once  k  is  known. 

Combining  Eq.  6  and  11  gives  a  relation  between  the  charging 
function  and  the  steady-state  hydrogen  permeation  current  density's, 29 

Equation  17  relates  the  charging  function,  4  exp(aTiF/Rr),  which  is 
solely  determined  by  the  conditions  at  the  charging  surface  of  the 
membrane  to  which  is  determined  by  the  transport  of  hydrogen 
through  the  membrane.  It  predicts  a  straight-line  relation  between 
the  charging  function  and  I’ooJ  from  the  slope  and  intercept,  one  can 
calculate  ^|  (=  zyZ'Cjj+),  and  k.  Then,  k  can  be  used  in  Eq.  1 1  to 
obtain  the  hydrogen  surface  coverage,  0fj,  as  well  as  k2  (Eq.  16).  The 
values  of  k  for  fire  different  membrane  thicknesses  can  be  used  to 
yield  ^des  described  above  (Eq.  15). 

Substiteting  for  from  Eq.  4  in  Eq.  9  and  rearranging,  one 
obtains^'^^ 


-  IS^L/D  [18] 

Equation  18  provides  for  the  calculation  of  the  hydrogen  concentra¬ 
tion  in  the  membrane  at  the  input  surface  once  k  and  0^  are  known, 
and  is  an  alternative  to  calculating  C®  from  Eq.  4  using  the  measured 
Zco  value. 


Experimental 

The  electrochemical  hydrogen  permeation  cell  was  similar  to  that 
used  by  Frumkin^®  and  by  Devanathan  and  Stachurski.^i  The  cell 
was  used  to  collect  data  on  both  the  HER  and  HAR  on  iron  mem¬ 
branes  of  different  membrane  thicknesses.  This  thickness  range  for 


Journal  of  The  Electrochemical  Society,  147  (8)  2959-2963  (2000) 
80013-4651(99)09-067-9  CCC:  S7.O0  ©  The  Elcctrochcinical  Society,  Inc. 


2961 


the  permeation  measurements  avoided  excessive  surface  control  and 
short-circuit,  e.g.,  grain  boundary,  permeation  of  hydrogen  while 
permitting  the  attainment  of  measurable  steady-state  permeation 
rates  which  were  achieved  in  less  than  1 0  min  for  all  runs.^®’^^  The 
iron  (steel)  membranes  with  thicknesses  of  0.10,  0.25,  0.50,  and 
0.85  mm  were  obtained  from  Goodfellow  with  the  following  com¬ 
position:  Mn,  0.3%;  Si,  0.1%;  C,  <0.08%;  S  <0.05%;  the  balance 
was  iron.^  They  were  polished  down  to  0.5  p.m  alumina  and  then  de¬ 
greased  in  acetone  and  washed  with  double-distilled  water.  They 
were  subsequently  annealed  in  pure  hydrogen  at  900°C  for  2  h  in  a 
tube  furnace  and  furnace-cooled  in  the  same  atmosphere.  The  exit 
surface  of  the  iron  membrane  was  coated  with  electrolessly  deposit¬ 
ed  palladium  using  Pallamerse  solution  before  the  membrane  was 
mounted  in  the  cell.  The  electrolyse  in  the  charging  (input)  compart¬ 
ment  of  the  cell  was  an  acidic  solution  of  0.1  N  H2SO4  +  0.9  N 
Na2S04  of  pH  1.8.  It  was  prepared  from  analytical  grade  chemicals 
and  double-distilled  water.  Before  admitting  the  solution  to  the  cell, 
it  was  pre-electrolyzed  at  3  mA  for  2  h  to  remove  impurities  that 
could  otherwise  affect  the  quality  of  the  data.  The  solution  was  sub¬ 
sequently  deaerated  with  hydrogen.  The  exit  compartment  of  the  cell 
contained  0.1  N  NaOH,  and  the  potential  of  the  Pd-coated  exit  sur¬ 
face  of  the  iron  membrane  was  set  at  0.150  mV(Hg/HgO)  using  a 
potentiostat  to  ensure  complete  oxidation  of  the  dissolved  hydrogen 
diffusing  through  the  membrane. 

To  minimize  the  production  of  hydrogen  traps  in  the  membrane 
dunng  the  measurements,  the  cathodic  current,  4,  at  the  charging 
surface  was  set  at  its  highest  value  and  then  decreased  stepwise  after 
reaching  within  minutes  the  steady-state  permeation  current  at  each 
value  of  Zj..  Similarly,  for  the  polarization  curve  measurements  in  the 
same  charging  compartment  of  the  cell,  the  cathodic  current,  z^.,  was 
decreased  stepwise  after  reaching  in  seconds  the  stationaiy  electrode 
potential  of  the  charging  surface  at  each  value  of  z^.  Further  details 
on  the  experimental  setup  are  reported  elsewhere.28,29,31,33 

Results  and  Discussion 

The  current-potential  relations  measured  on  the  charging  side  of 
each  membrane  gave  a  Tafel  line  with  a  slope  of  120  mV,  indepen¬ 
dent  of  the  membrane  thickness.  This  corresponds  to  a  value  of  a  = 
0.5.  Figure  2  shows  plots  of  Zoo  v^.  for  membranes  of  different 
micknesses.  These  data  plot  well  as  straight  lines  passing  through 
the  point  of  origin,  in  accordance  with  Eq.  16,  indicating  that  the 
mechanism  of  the  HER  is  coupled  Volmer  discharge-Tafel  recombi¬ 
nation.  The  slopes  of  the  lines  in  Fig.  2  decrease  with  increasing 
membrane  thickness,  as  predicted  by  Eq.  12  and  16. 

Fi^re  3  shows  the  straight-line  relations  obtained  between  the 
charging  function,  z^  exp(Q''qF7R7^,  and  the  steady-state  hydrogen 
permeation  current  density,  z«,  for  different  membrane  thicknesses. 
Table  I  lists  the  values  of  z^  estimated  from  the  intercept  of  these 
plots  m  accordance  with  Eq.  17,  the  slopes  of  the  lines  in  Fig.  3,  and 
the  kinetic-diffiision  constant,  k,  obtained  from  these  slopes  (Eq.  17) 
on  iron  membranes  of  different  thicknesses  in  0. 1  N  H2SO4  +  0,9  N 
Na2S04.  Table  I  shows  that  as  the  membrane  thickness  increases,  the 
value  of  the  kinetic-diffusion  constant,  k,  decreases.  It  further  shows 


Figure  2.  The  relation  between  the  steady-state  hydrogen  permeation  current 
density,  and  square  root  of  the  recombination  current  density,  obtained 
on  iron  membranes  of  different  thicknesses  in  0. 1  N  HjSO.  -I-  0  9  N  Na,SO. 
(pH  1.8).  ^  •  2  4 


1.60 

1.50 

1.40 

1.30 

1.20 

1.10 
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0.90 

0.80 


L,  pA  cm'* 

Figure  3.  The  relation  between  the  charging  fiinction,  ex^iar\FlRr),  and 
hydrogen  permeation  current  density,  z„,  for  iron  membranes 
of  different  thickness  in  O.I  N  H2SO4  +  0.9  N  Na2S04. 

that  the  exchange  current  density  does  not  change  with  the  mem¬ 
brane  thickness. 

The  reciprocal  of  the  values  of  the  kinetic-diffusion  constant,  k, 
was  plotted  in  Fig.  4  vs.  the  membrane  thickness,  L,  to  evaluate  Ic  u’ 
and  according  to  Eq.  15.  Figure  4  shows  a  straight-line  relation 
with  a  positive  slope  and  positive  intercept  with  they  axis,  in  accor¬ 
dance  with  Eq.  15  and  the  IPZ  analysis.  The  diffusivity  of  hydrogen 
within  the  membranes  was  obtained  using  the  breakthrough  time.^^ 
A  value  ofD  =  1.2  ±  0.02  X  10“^  cm^  s~^  was  found  for  all  thick¬ 
nesses.  This  value  cornpares  well  with  reported  values  for  the  diffu¬ 
sivity  of  hydrogen  in  iron.^'^  It  also  reveals  that  the  diffusivity  was 
not  a  function  of  the  membrane  thickness.^^ 

From  the  slope  and  intercept  of  Fig.  4,  we  calculated  values  of 
4bs  =  2.0  X  10  mol  cm  ^  s~^  and  =  1.9  X  10“^  cm  s“^ 
This  value  of  can  be  compared  to  the  values  of  1-3  X  lO^^  cm 
s  measured  in  solutions  of  pH  13  on  high  purity  iron  samples 
ranging  from  50  to  900  jxm  thick,^^  1 0*“2  cm  s“  ^  measured  in  H2SO4 
on  Armcoffon  ranging  in  thickness  from  0.026  to  0.129  cm, 23  and 
10  cm  s  ^  measured  on  commercially  pure  iron  (99.7%)  in  0  2  N 
NaOH.36  The  value  of  =  2.0  X  lO"^® 
compared  to  the  values  reported  recently  by  Zhang  et  al  in  the 
fo^  of  the  product  Qk^^^  ranging  from  10“^^  to  10"^^  mol  cm-^ 
s  .  There  are  also  the  values  reported  by  Kim  and  Wilde^^  and  Sub- 
ramanyan.  Both  of  these  works  used  absolute  values  of  coverage 
and  hence  obtained  units  of  s"^  namely,  3-5  and  1  5  s“^ 

respectively.  To  convert  these  values  of  from  units  of  s'^  intj 
units  of  mol  cm  ^  $  \  one  must  multiply  them  by  the  monolayer 
capacity  of  hydrogen  adsorbed  on  iron.  A  value  of  2  X  10'^  mol 
cm  was  used  by  Kim  and  Wilde,^®  which  is  essentially  an  arbitrary 
value.  Using  this  value,  their  rate  constants  of  3-5  and  1  5  s"^  give 
values  of  6-10  X  IQ-^  and  3  X  10“^^  mol  cm“2  s“i.  Using  = 
2X10  mol  cm  s  ^  and  -  1.9  X  10”^  cm  s“^,  the  contri¬ 
bution  of  surface  control  can  be  estimated  from  Eq.  15.  The  magni¬ 
tudes  of  the  first  term  are  approximately  7  and  38%  of  the  sum  of 
both  terms  on  the  right  side  of  Eq.  15  for  the  0.85  and  0,10  mm 
membranes,  respectively  (Table  I). 

The  hydrogen  surface  coverage,  63,  was  calculated  from  Eq.  11 
using  the  k  values  obtained  for  the  membranes  of  different  thick¬ 
nesses.  Figure  5  shows  that  the  magnitude  of  Oh  increases  with  more 


Table  I.  Values  of  -i'^/Fk,  and  k  (Eq.  17)  and  the  %  surface 
control  for  Iron  membranes  of  different  thicknesses  in  0.1  N 
HjS04  +  0,9  N  NajS04  (pH  1.8). 


L  (mm)  (A  cm"^  _ ^ 

Fk 


/t(molcm-2s-')  %  surface 


control 


0.10 

2.0  X  10“® 

0.32 

6.6  X  10"'> 

0.25 

2.0  X  10“® 

0.40 

5.4  X  10-“ 

0.50 

2.0  X  10“^ 

1.20 

2.1  X  10-“ 

0.85 

2.0  X  ]0~® 

1.45 

1.4  X  10-“ 

38 

20 

10 

7 
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Figure  4.  The  relation  between  the  reciprocal  of  the  kinetic-diffiision  con¬ 
stant,  k,  and  the  membrane  thicknesses,  L,  for  iron  in  0.1  N  H2SO4  +  0.9  N 
Na2S04. 


Figure  5.  The  relation  between  the  hydrogen  surface  coverage,  0^,  and  the 
electrode  potential,  E,  obtained  on  iron  membranes  of  different  thicknesses 
in  0.1  N  H2SO4  +  0.9  N  Na2S04,  and  the  corresponding  hydrogen  concen¬ 
tration,  C°,  in  the  iron  membrane  at  the  input  surface. 


negative  E  values,  in  accordance  with  its  expected  dependence  on 
the  overpotential  This  trend  is  well  documented  in  the  litera- 
ture.^^’'*^  It  follows  that  0jj  is  also  independent  of  the  membrane 
thickness,  as  indicated  by  the  data  in  Fig.  5.  These  values  of  0^^  were 
obtained  from  the  model  assuming  the  Langmuir  adsorption  iso¬ 
therm.  Others  have  obtained  similar  values  for  0^  on  iron  and  steels 
also  using  the  Langmuir  adsorption  isotherm.'^o,46  Support  for  the 
use  of  the  Langmuir  adsorption  isotherm  also  comes  from  the  results 
using  the  Fmmkin  isotherm,  which  show  that  the  interaction  para¬ 
meter  is  negligible  (of  the  order  of  0.2).'**^  From  the  values  of  the 
hydrogen  surface  coverage,  the  hydrogen  concentration  in  the  iron  at 
the  input  surface,  C°,  can  be  obtained  from  Eq.  18.  These  values  are 
also  shown  in  Fig.  5  and  are  of  the  same  order  of  magnitude  as 
values  obtained  from  the  diffusion  law,  Eq.  4  (using  the  data  in 
Table  II  and  the  measured  D  values). 

Conclusions 

Applying  the  IPZ  analysis  to  the  measured  steady-state  hydrogen 
permeation  data  for  different  values  of  charging  current  and  mem¬ 
brane  thickness,  the  specific  rate  constants  of  hydrogen  absorption 
and  desorption  were  obtained  for  iron  in  an  acid  sulfate  solution.  The 
analysis  utilized  a  corrected  Eq.  15  from  that  given  in  the  original 
work.28,29  Using  these  results,  the  amount  of  surface  control  during 
the  hydrogen  absorption-permeation  process  could  be  evaluated  for 


each  membrane  thickness.  The  analysis  yielded  also,  as  previously 
obtained  for  a  single  membrane  thickness,  the  degree  of  hydrogen 
surface  coverage,  Gy,  exchange  current  density,  and  the  rate  con¬ 
stants  of  the  HER.  It  also  yielded  the  important  kinetic- diffusion  con¬ 
stant  k  derived  previously^"^  which  is  proportional  to  the  k'  constant 
(k  =  k'D/L)  presented  previously, ^8.29  which  combines  the  rate  con¬ 
stants  of  absorption  and  desorption  at  the  charging  surface  and  the 
transport  properties  of  the  membrane.  For  the  iron/acidic  sulfate  sys¬ 
tem,  the  IPZ  analysis  gives  =  2  X  10”^®  mol  cm“^  s~^  and 
^des  =  L9  X  10“^  cm  s"^,  in  fair  agreement  with  the  few  values 
obtained  by  other  methods  in  the  literature.  The  hydrogen  surface 
coverage  was  a  function  of  the  applied  cathodic  potential  in  accor¬ 
dance  with  the  literature,  but  not  of  the  sample  thickness  (Fig.  5). 
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Cjj+  hydrogen  ion  concentration,  mol  cm^^ 
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ki  discharge  rate  constant,  cm  8“^ 

}c2  recombination  rate  constant,  mol  cm“^  s“  * 

/Tabs  absorption  rate  constant,  mol  cm“^  s”  ^ 

^des  desorption  rate  constant,  cm 

k  kinetic-diffiision  constant,  mol  cm"^  s“* 

R  gas  constant,  8.314  J  mol“^ 

T  absolute  temperature,  K 
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a  transfer  coefficient,  dimensionless 

Tj  cathodic  overvoltage,  V 

hydrogen  surface  coverage  at  equilibrium,  dimensionless 
6h  hydrogen  surface  coverage,  dimensionless 
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An  Analysis  Procedure  for  Hydrogen  Absorption  under 
Frumkin  Adsorption  Conditions 
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This  paper  communicates  a  more  genera!  derivation  of  ihc  lycr-Pickcring-Zamcnzadch  (IPZ)  analysis  of  hydrogen  absorption.  It 
is  applicable  to  mctal/clcctrolytc  .sy.stcms  under  Frumkin  adsorption  conditions  for  the  discharge-recombination  process  oHiydro- 
gen  evolution,  as  well  as  to  Langmuir  adsorption  conditions  covered  by  the  original  derivation.  As  such,  it  does  not  require  prior 
knowledge  of  the  dimcnsionics.s  factor/ that  describes  the  deviation  from  ideal  Langmuir  behavior.  Titis  analysis  appropriately 
reduces  to  the  original  fPZ  analy.sis  for  Langrrniir  ad.sorption  conditions  for  which  /  =  0.  It  yields  ihc  value  of  (/')  as  well  as  the 
other  unknowns  in  the  system:  the  discharge  (k,)  and  recombination  (A^)  rate  constants,  the  exchange  current  density  (/p),  the 
kinetic-diffusion  constant  (A),  the  hydrogen  surface  coverage  (fin)-  hydrogen  concentration  in  the  metal  at  the  charging 

surface  (C®). 
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The  Frumkin  adsorf'>tion  isotherm  assumes  that  the  free  energy  of 
ad.sorption  of  a  specic.s  decreases  with  coverage  according  to  the 
following  equation* 

AG°  ==  + /RTO  [1] 

where /is  a  dimensionless  factor  that  de.scribes  the  deviation  from 
the  ideal  Langmuir  behavior.  The  fRT  quantity  is  sometimes  called 
the  rate  of  change  of  the  standard  free  energy  of  adsorption,  AGJJ 
and  AGl  are  the  standard  free  energy  of  adsorption  at  a  certain 
coverage,  0,  and  at  zero  coverage,  0  =  0,  respectively,  R  is  the  gas 
constant,  and  T  is  the  absolute  temperature. 

The  Langmuir  adsoiption  isothemi  assumes  the  .standard  free  en¬ 
ergy  of  adsorption  is  independent  of  coverage  (/  =0).  This  iso¬ 
therm  is  usually  applicable  at  low  values  of  coverage  (0  <  0.1)  and 
at  coverages  close  to  unity  (0  >  0.9).'^’^  In  some  instances  it  has 
been  found  to  be  applicable  also  for  mtcrmediaie  coverages  of  hy¬ 
drogen  on  iron  surfaces.*"'*^ 

The  Tycr-Pickering-Zamanzadeh  (IPZ)  analysis’’^  developed 
three  major  rclationslups  between  the  hydrogen  evolution  reaction 
(HER)  and  the  hydrogen  ad.sorption  reaction  (HAR).  The  first  one  i.s 
between  the  .steady-state  hydrogen  permeation  current  density,  U , 
and  the  hydrogen  recombination  current  density,  if.  This  relation  is 
shown  in  Eq.  2.  Another  rclalioaship  is  between  tlic  charging  fiinc- 
tion,  /<.  e-xpCcictH).  and  /*  which  is  .shown  in  Eq.  3.  The  last  one  i.s 
between  and  the  hydrogen  coverage,  0h,  shown  in  Eq.  4 


I'c  exptr/aT))  =  1  ~  W 

ire  =  Fkdn  W 

where  is  the  charging  current  density.  /'  =  F/ciCh-j.  = 

(1  ~  Gf,),  /’o  is  the  exchange  current  density  of  the  HER  and  Of,  is 
the  hydrogen  .surface  coverage  at  equilibrium.  C,‘,  is  the  hydrogen 
ion  conccnU'ation  in  the  electrolyte,  a  =  F/RT  where  F.  R.  and  T 
have  their  usual  meanings,  a  is  the  transfer  coefficient,  and  h  is  the 
overpotential  of  the  HER.  k  is  the  kineiic-dilfusion  con.stant  and  is 
defined  as 
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k  =  — ^  [5] 

where  is  the  rate  con.stant  for  hydrogen  absorption  into  the  metal 
at  the  charging  surface,  is  the  desorption  rate  constant  for  the 
hydrogen  moving  in  the  opposite  direction,  i.e.,  from  the  absorbed  to 
ad.sorbcd  state,  D  is  the  IiycU'ogen  diffusion  coefficient  within  the 
membrane,  and  L  i.s  the  membrane  thickness. 

Equation  5  is  a  modified  version  of  Eq.  4a  in  tlie  original  IPZ 
anuIysLs^  and  was  recently  derived  by  Rama.subramanian  c/r/A  *** 
More  details  about  the  IPZ  analysis  can  be  found  elsewhere.^'®  For 
the  local  equilibrium  of  the  intermediate  absoq^tion-desoq^tion  reac¬ 
tion,  the  concentration  of  hydrogen  inside  the  metal  is  obtained  as^ 

C®  =  A'OhL/O  [6] 

Experimental 

The  electrochemical  hydrogen  permeation  cell  vva.s  similar  to  that 
u.sed  by  Devanathan  and  Stachurski,**  The  cell  was  used  to  collect 
data  on  the  HER  and  HAR  on  iron  membranes  of  thickness  0.25  mm 
in  acidic  solution  of  I  M  Na2S04.  acidified  to  pH  2  using  H2SO4 
(blank  solution)  and  blank  solution  with  a  small  addition  of  thio.sul- 
fatc  (S2O3").  Thiosulfate  is  known  to  be  a  promoter  for  hydrogen 
absorption  by  iron  which  docs  not  satisfy  the  requirement  of  the 
original  IPZ  analysis. The  membrane  was  coated  with  palladium 
only  on  its  exit  side.  Further  details  about  the  experimental  setup  can 
be  found  elsewhere;^  • 

Analysis  Development 

Using  the  Volmer  dischargc-Tafcl  hydrogen  adatom  combination 
process  of  hydrogen  evolution  at  appreciable  oveipotentials  (so 
there  is  negligible  hydrogen  oxidation),  whicii  was  the  basis  for  liie 
original  IPZ  analysis,^'^  the  HER  proceeds  as  follows 


*1 


+  M  -F  e"' 

[7] 

,  +  H;  +  2M 

[8] 

Some  of  the  ad.sorbcd  hydrogen  atoms,  will  absorb  into  the 

iron,  and  quickly  cstabli.sh  the  back  equilibria 

^  +  M  [9] 

The  absorbed  liydrogcn  diffuse.^  through  the  iron  membrane. 
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Tabic  1.  Summary  nf  major  relations'  developed  by  the  original  IIY,  aiirilysi.s  and  Hr*  generalized  (Frunikin)  IVZ  analysis  derived  in  tliis  paper. 
Frumkin  adsorption  isoihcnn  Langmuir  adsorption  Isotherm 

ACTS  =  Ar;g  +  /RTn  [I]  ic;  =  ac; 

Generalized  IPZ  analysis  Origimil  IPZ  analysis 


[7-1 

f.l] 


In  writing  the  rale  equations  in  view  of  the  Faimkin  adsorption 
conditions,  one  has  for  I  he  rate  of  proton  discliargc  and  hydrogen 
evolution^'^'*'^ 


cxp(i/aT]) 


[19] 


i\  =  “  0„)  oxp(-aaT\)  exp(-a/0//)  [10] 

t'r  =  ‘^>^.p(2a/0//)  [II] 

where  and  arc  the  current  densitie.s  of  proton  discharge  (Fq.  7) 
and  hydrogen  adatom  combination  (Eq.  8),  respectively. 

Inserting  0}j  from  Eq.  4  into  Eq.  11 


Taking  the  square  root  of  Eq.  12 

V?-T' 


Dividing  both  sides  ofEq.  13  by  and  taking  the  natural  logarithm 
(In) 


This  relation,  which  wa.s  also  denved  by  Iyer  ef  ai^  in  a  similar 
form  but  with  different  slope  and  intercept,  can  be  used  for  deter¬ 
mining  the  rate  constants  when  the  value  of/ is  known  or  can  oth¬ 
erwise  be  cstimated.*’^^'*'*  An  independent  evaluation  of/  however, 
is  not  an  easy  task.  Tn  the  following  derivation,  we  arrive  at  anotlicr 
relationship  which,  when  applied  to  the  steady-state  permeation 
data,  can  yield  the  rate  constants  without  a  prior  knowledge  of  the 
value  of/  In  fact,  the  value  of/ is  a  result  of  the  analysis. 

Rcan’anging  Hq.  10  and  solving  for  exp(a/flu) 


expCa/On)  = 


/o(l  -  Bh) 
exp(£/ot'q) 


D5] 


Taking  the  square  root  of  Eq,  1 1 

\fFk^  •  Oh  •  cxpCo/Oh) 


Substituting  Eq.  15  into  Eq.  16 


^ff,  =  v'«3  •  Oh  • 


'od  -  Oh) 

/‘c  e.xpCocxT)) 


[16] 


[17] 


Rcarraitging  and  inserting  for  0J^  from  Eq.  4  into  Eq,  17 


cxpiaixTi)  =  •  — 


[IS] 


Equations  14  and  19  can  be  simplified  to  Eq.  2  and  3  if  the  Lang¬ 
muir  adso>i>lion  conditions  are  assumed  instead  of  the  Frumkin  ad¬ 
sorption  conditions.  In  other  words,  the  derived  IPZ  analysis  in  this 
paper  reduces  to  the  original  Langmuir  IPZ  analysis’  if  the  standard 
free  energy  of  adsorption  is  assumed  to  be  independent  of  the  hy¬ 
drogen  coverage,  /.e.,  /  “  0.  Thus,  the  original  IPZ  analysis  is  a 
special  case  of  the  IPZ  analysis  in  this  paper,  in  the  same  way  the 
Langmuir  isothenn  is  a  special  case  of  the  Frumkin  isotherm.  These 
relations  are  summarized  in  Table  L  ^ 

Plots  of  both  ln(\////x.)  and  v.y.  ty,  should 

give  straight  lines  according  to  Eq.  14  and  19,  respectively. 

The  procedure  of  calculating  the  unknowns  in  the  system  is 

1.  From  the  slope  and  intercept  of  Eq.  19,  k  can  be  calculated. 

2.  From  the  slope  of  Eq.  14  and  knowing  k  from  step  I,/ can  be 
calculated,  after  using  the  charging  (Tafel)  daui  to  obtain  the  value 
of  a. 

3.  From  the  intercept  of  Eq.  14  and  knowing  k  from  step  1, 
can  be  calculated. 

4.  From  the  intercept  of  Eq.  19  and  know'ing  k  from  step  1  and  /c, 
from  step  3,  can  be  calculated. 

5.  From  Eq.  4  and  knowing  /c  from  step  1,  the  hydrogen  surface 
coverage  (6).|)  can  be  calculated. 

Results  and  Discussion 

Figure  1  shows  the  relationship  between  the  steady-state  hydro¬ 
gen  permeation  current  density,  /* ,  and  the  square  root  of  the  hy¬ 
drogen  recombination  current  density,  v'/p,  whore  if  is  obtained  as 
the  difference  between  the  impressed  i^  value  and  the  measured 
value.’  The  figure  show'.s  a  straight  line  passing  by  the  origin  for  the 
blank  (/.e.,  without  S2O3")  .solution.  This  agrees  with  Eq.  2  and 
indicates  that  the  HER  occurs  according  to  the  reaction  .sequence  in 


Fijjurc  L  The  relation  between  the  measured  .steady-state  hydrogen  perme¬ 
ation  current.  and  the  square  root  of  the  recombination  current,  v^- 
Blank  =  I  M  Na2S04.  pH  2.  TS  =  S^Oj"'. 


Dividing  both  side.s  of  Eq.  18  by  /. 
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Figure  2.  The  rclution  between  tijc  charging  funciion,  i\.  expiaa-r]).  and  ihc  3-  The  v;iriation  of  the  hydrogen  surface  coverage,  0,^.  with  the 

measured  steady-state  hydrogen  permeation  current,  4,  for  the  blank  (i  M  electrode  polcntsal  for  the  blank  (1  M  Na2S04,  pH  2.)  solution. 

Na2S04.  pH  2J  solution. 


Eq.  7  and  8,  that  Eq.  7  is  not  near  equilibrium,  hydrogen  oxi¬ 
dation  is  negligible,  and  that  the  IPZ  analysis  can  be  applied  to  llic 
steady-state  hydrogen  penneation  data7‘^ 

However,  in  the  presence  of  thiosulfate,  the  line,  although  linear 
CFig.  1),  docs  not  pass  by  the  origin  and  has  a  finite  intercept  with 
the  y  axis.  Therefore,  the  original  IPZ  analysis  cannot  be  used  to 
evaluate  the  steady-state  hydrogen  permeation  data  in  the  presence 
of  thiosulfate,  as  di.scussed  elsewhere*  based  on  the  analysis  of  hy¬ 
drogen  permeation  data  for  H^S  additions  to  acidic  chlorate  solu¬ 
tions.  We  therefore  used  these  results  as  an  example  to  lest  the 
applicability  of  the  analysis  which  is  based  on  the  Frunikin  ap¬ 
proach,  in  spite  of  the  fact  that  there  is  no  independent  evidence  that 
the  Frumkin  approach  Is  a  better  description  than  the  Langmuir  ap¬ 
proach  for  H  ad.sorption  under  these  conditions.  Furthermore,  we 
ignored  Any  effect  that  the  sulfur  species  may  have  had  on  the  den¬ 
sity  of  available  hydrogen  adsorption  sites  on  the  .surface,  for  (he 
sake  of  an  approximation. 

Figure  2  shows  the  relationship  berw’een  the  charging  function 
[/‘e  exp(r/a7i)]  and  the  measured  steady-state  hydrogen  permeation 
current  demsiiy,  for  the  blank  solution.  The  figure  shows  a 
straight  line  with  a  negative  slope  and  positive  intercept  in  accord 
with  Eq.  3.  The  rate  constants  of  the  HER,  and  k2,  the  kinetic- 
diffusion  constant  of  the  HAR.  k\  the  exchange  current  density,  /^, 
and  the  corresponding  hydrogen  concentration  just  inside  the  metal 
at  the  charging  surface,  C“,  were  calculated  for  thc.se  data  obtained 
in  the  blank  (thiosulfate-free)  solution  using  tlie  original  IPZ  analy¬ 
sis.  Results  arc  summarized  in  Table  11.  Hydrogen  coverages.  On, 
and  hydrogen  concentrations,  C\  at  different  electrode  potentials 
(£■)  were  calculated  and  the  results  arc  show'ii  in  Fig.  3  and  4,  re¬ 
spectively,  Details  on  these  calculations  can  be  found  elsewhere.’*^ 
The  steady-state  hydrogen  permeation  data  for  both  the  blank 
solution  and  thiosulfate  solution  were  also  tested  using  the  general¬ 
ized  IPZ  analysis  derived  in  this  paper.  Figure  5  shows  a  straight 


Table  II.  Value.s  of  tilfferent  constants  obtained  by  applying  the 
original  (Langmuir)  IPZ  and  the  generalized  (Frumkin)  IPZ 
analyses  to  measured  steady-state  permeation  data  for  iron  in  the 
blank  solution  (I  M  Na2.S04.  pH  2)  or  in  the  blank  +  I  mM 
solution. 


Generalized  1P2 

Original  IPZ  Generalized  IPZ  (blank  +  I  inM 
Variable  (blank)  (blank)  S2O5) 


k,  mol  cjir^  s”' 
^2,  mol  cm*'^  s" ' 
,  p.A  cm~“ 

/tj  ,  cm  s'* 

f 


1.02  X  10““’  1.02  X  lO*"'” 

5.02  X  4.63  X  JO"*' 

1.17  1.17 


1.21  X  10"*'  1.21  X  10'^’ 

0  (assumed)  0.20 


2.68  X  10“^'“ 
2.28  X  lO’’** 
34.4 

3.56  X  )0‘  ^ 
9.90 


line  relation sliip  between  and  the  measured  steady-state 

hydrogen  permeation  current  density,  in  accord  with  Eq.  14. 
Figure  6  shows  a  linear  relation  between  ( and 
the  measured  stciidy-state  hydrogen  permeation  current  density,  /y,, 
in  accord  with  Eq.  19.  The  procedure  summarized  here  was  used  to 
calculate  the  unknowns.  The  Tafel  slope  was  120  mV,  which  corre¬ 
sponds  to  a  value  of  ot  =  0.5.  The  results  arc  shown  in  Table  U.  The 
calculated  hydrogen  coverage  using  the  generalized  IPZ  analyses  is 
also  shown  in  Fig.  3.  Since  it  is  unknown  to  what  extent  the  .sulfur 
species  in  the  solution  may  have  affected  the  density  of  available 
adsorption  sites  for  hydrogen,  the  values  for  the  thiosulfate  solution 
in  Tabic  11  could  be  limiting  ones. 

Results  using  the  two  different  IPZ  analyses  were  found  to  be 
similar  when  the  requirements  of  both  the  original  and  generalized 
IPZ  analyses  were  fulfilled.  Only  the  generalized  IPZ  analysis  can 
be  applied  when  the  requirement  of  the  original  iPZ  analysis  is  not 
fulfilled  (/>.,  the  relation  between  v.v.  v7^  is  not  linear  and/or  docs 
nor  pass  by  the  origin,  e.g.,  the  blank  +  I  mM  TS  line  in  Fig.  1). 

The  factor/ was  assumed  to  be  zero  in  the  derivation  of  the 
original  iPZ  analysis.  The  generalized  IPZ  analysis  yields  a  value 
close  to  zero  (./'  =  0.20)  for  the  blank  solution,  whereas  it  yields  a 
value  of/  =  9.9  for  the  blank  solution  containing  I  mM  thiosulfate 
solution.  The  good  correlation  between  the  original  and  generalized 
fPZ  analyses  (Fig.  3)  along  with  the  low  (essentially  zero)  calculated 
value  of/ support  the  existence  of  the  Langmuir  condition  in  the 
case  of  the  blank  solution. 


Conclusions 

The  original  IPZ  analysis  has  been  made  more  gcncrdl  so  as  to  be 
applicable  to  steady-state  hydrogen  penneation  data  for  a  wide 
range  of  metal/el ectrolylc  systems  undergoing  the  HER.  This  new 
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Figure  4.  The  variation  of  the  hydrogen  concentration  just  inside  the  cliarg- 
ing  surface,  0\  with  the  electrode  potential  for  the  blank  (I  M  Na^SO^.  pH 
2)  solution. 
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Figure  5.  The  relation  between  steady-state  permeation 

current  density,  iy,.  Blank  =  I  M  Na2S04.  pM  2,  TS  =  S^Oj  . 


IPZ  analysis  improves  upon  an  earlier  version  derived  for  Frunikin 
adsorption  conditions.  It  is  more  u-sefuJ  than  the  earlier  version  be¬ 
cause  it  does  not  require  a  prior  knowledge  of/  as  its  predecessor 
did.  Rather,  the  IPZ  analysis  derived  in  this  paper  yields  the  value  of 
/  for  those  metal/eleclrolyte  systems  that  meet  the  conditions  for 
application  of  the  TPZ  analysis  given  clscvvlierc.*^'^  in  addition,  it 
yields  the  rate  constants  and  ^2)  exchange  current  den¬ 

sity,  to,  of  the  HER,  the  kinetic-diffusion  constant,  k,  of  the  HAR, 
hydrogen  surface  coverage,  and  the  hydrogen  concentration  just 
inside  tlie  metal  surface,  C®,  just  as  the  original  IPZ  analysis  does 
for  Langmuir  adsorption  conditions. 

The  generalized  (Fnimkin)  IPZ  analysis  can  be  simplified  to  the 
original  (Langmuir)  IPZ  analysis  when  the  standard  free  energy  of 
adsorption  is  independent  of  the  surface  coverage,  =  0.  The 
results  obtained  for  a  system  exhibiting  the  Langmuir  adsorption 
conditions  are  tlie  same  using  either  the  original  or  the  generalized 
IPZ  analyses.  However,  the  generalized  IPZ  analysis  must  be  used 
for  the  systems  which  do  not  meet  the  Langmuir  conditions,  i.e.,  for 
systems  in  which  the  measured  steady-state  hydrogen  permeation 
data  do  not  satisfy  Eq.  2  and/or  Eq.  3,  e.g.,  see  Fig.  4b  and  6b  of 
Ref.  8  for  the  HER  in  a  pH  2  chlorate  solution  with  small  addition  of 
H2S.  Still,  the  applicability  of  these  analyses  in  more  complicated 


£»  (iA  cm*’ 

Figure  (u  The  relation  between  ( exptrta-n)  and  the  steady-state 
permeation  current  den.slty,  .  Blank  ==  1  M  Na2S04.  pH  2.  TS 
=  S,Oj  \ 


electrolytes,  in  particular  when  the  availability  of  hydrogen  ad.sorp- 
tion  sites  could  be  affected  by  other  solution  species,  needs  to  be 
taken  into  consideration. 
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L  membntne  thickness,  cm 
R  gas  constant,  8.314  J  (mo!  K)"' 

T  absolute  tempemture,  K 

Greek 

a  transfer  cocfHcicnt,  dimensionless 
71  oveqiotcniial  of  the  HER,  V 
f)j.|  hydrogen  surface  coverage,  dimensionless 
Oh  hydrogen  surface  coverage  ;U  equilibrium,  diincnsi<inless 
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Effect  of  Thiosulfate  and  Sulfite 

on  the  Permeation  Rate  of  Hydrogen  Through  Iron 

M.H.  Abd Elhamid, *  B.G,  Ateya,  **  K.G.  Well,  and  H.W.  Pickerlng'^^*^ 


ABSTRACT 


BJfects  of  thiosulfate  (S20ii  and  sulfite  (SOt)  on  hydrogen 
evolution  and  permeation  within  iron  have  been  studied  in 
acid  and  neutral  solutions.  Results  show  that  thiosulfate 
increases  the  hydrogen  permeation  rate  and  depolarizes  the 
hydrogen  evolution  reaction  at  concentrations  as  low  as 
0. 1  mM,  in  acid  or  neutral  media.  Sulfite  showed  similar  and 
somewhat  less  pronounced  effects  in  the  acid  medium  and 
no  effect  on  either  process  in  the  neutral  medium.  Results  are 
discussed  in  light  of  the  stability  of  the  various  sulfur  species 
and  their  effects  on  the  rates  of  both  processes. 

KEY  WORDS:  environmentally  assisted  corrosion,  hydrogen 
evolution,  hydrogen  overvoltage,  hydrogen  permeation,  steel 
sulfide  stress  corrosion,  thiosulfate 

INTRODUCTION 


The  degrading  effects  of  sulfur-containing  (thio)  com¬ 
pounds  on  the  integrity  of  metals  and  alloys  have 
been  well  recognized.  Of  these  compounds,  thiosul- 
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fate  is  commonly  encountered  In  the  pulp,  paper,  and 
oil  and  gas  industries  and  as  a  product  of  sulfate- 
.  reducing  bacteria.  For  this  reason,  extensive  studies 
have  been,  conducted  by  various  authors  to  evaluate 
and  understand  the  effects  of  thiosulfate  on  crevice 
corrosion,  pitting  corroslon,^*^  stress  corrosion 
cracking,®'®  and  anodic  dissolution  of  metals  and 
alloys.*^'®  Furthermore,  some  authors  have  considered 
the  electrochemistry  of  thiosulfate®  and  its  adsorp¬ 
tion  behavior. 

Generally,  it  is  accepted  that  thiosulfate  has  a 
strong  promoting  effect  on  the  anodic  dissolution  of 
metals.'^’®  However,  measurements  of  its  effects  on 
hydrogen  absorption  have  been  less  systematic  and 
less  conclusive. 

The  objective  of  this  paper  was  to  evaluate  the 
effects  of  thiosulfate  and  sulfite  on  the  permeability 
of  electrolytic  hydrogen  within  iron.  Sulfites  were  In¬ 
cluded  in  this  study  since  thiosulfate  decomposes 
readily  to  give  sulfite  species  and  colloidal  sulfur  in 
mild  acidic  media  similar  to  those  media  encountered 
in  pits,  crevices,  and  cracks. 

EXPERIMENTAL  PROCEDURES 


An  electrochemical  hydrogen  permeation  cell^^*^’ 
was  used  In  this  study  to  collect  data  on  both  hydro¬ 
gen  evolution  and  permeation  within  iron  (steel) 
membranes  in  sulfate  solutions  containing  different 
thiosulfate  or  sulfite  concentrations.  Figure  1  shows 
the  position  of  the  iron  membrane  as  It  was  mounted 
between  the  two  parts  of  the  cell.  Steel  membranes 
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with  0,25  mm  thicknesses  were  obtained  from  a 
commercial  source  with  the  following  composition: 
0.3%  Mn,  0.1%  Si,  <  0.08%  C,  <  0.05%  S,  and  bal. 
Fe.  They  were  polished  successively  to  0.05-pm 
alumina  (AI2O3),  degreased  in  acetone  (CH3COCH5). 
washed  with  double-distilled  water,  and  subse¬ 
quently  annealed  in  pure  hydrogen  at  900°C  for  2  h 
in  a  tube  furnace,  followed  by  a  furnace  cool  in  the 
same  atmosphere.  Solutions  of  0.05  M  sulfuric  acid 
(H2SO4)  +  0.45  M  sodium  sulfate  (Na2S04,  pH  1.8) 
and  0.5  M  Na2S04  (pH  7)  were  used  throughout  this 
study.  All  solutions  were  prepared  from  analytical 
grade  chemicals  and  double -distilled  water.  Before 
admitting  the  solutions  to  the  cell,  they  were  pre- 
electrolyzed  using  platinum  electrodes  for  2  h  under 
deaeration  using  hydrogen.  These  were  then  used  as 
blanlc  solutions  to  which  different  concentrations  of 
thiosulfate  or  sulfite  were  added.  The  solutions  were 
subsequently  deaerated  with  hydrogen.  One  surface 
of  the  membrane  was  coated  with  palladium  before 
the  membrane  was  mounted  in  the  cell.  Areas  of  the 
input  (charging)  and  exit  (palladium-coated)  surfaces 
of  the  membranes  were  0.8  cm^.  The  palladium- 
coated  surface  was  held  at  0.150  mV  [mercuiy/mer- 
curic  oxide  [Hg/HgO])  in  0.05  M  sodium  hydroxide 
(NaOH)  to  ensure  complete  oxidation  of  the  atomic 
hydrogen  diffusing  through  the  membrane.  The  other 
(charging)  surface  of  the  membrane  faced  the  thio¬ 
sulfate  contaminated  solution  and  was  subjected  to 
hydrogen  charging  under  cathodic  polarization.  Fur¬ 
ther  details  about  the  experimental  arrangement  are 
reported  elsewhere. 

RESULTS  AND  DISCUSSION 


The  cathodic  charging  current  is  supported  by 
the  hydrogen  evolution  reaction  (HER)  taking  place  at 
the  input  surface  of  the  membrane,  l.e.: 

H*(aq)  +  e--^Had,  .  (D 

•  (2) 

where  H^ds  refers  to  a  hydrogen  atom  adsorbed  on  the 
metal  surface.  Reactions  (1)  and  (2)  represent  the 
discharge-recombination  mechanism  of  the  HER.  A 
fraction  of  this  adsorbed  hydrogen  (H^aJ  gets  ab¬ 
sorbed  into  the  lattice  of  the  metal  (Habs)  through  a 
surface  reaction: 


^ads  ^  ^abs  (3) 

Reaction  (3)  (hydrogen  absorption  reaction  [HAR]) 
leads  to  the  establishment  of  a  finite  concentration 
(C®)  of  hydrogen  atoms  absorbed  just  beneath  the 
charging  surface  of  the  membrane.  It  is  this  ab¬ 
sorbed  hydrogen  that  diffuses  through  the  membrane 
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FIGURE  1.A  schematic  illustration  showing  the  conditions  at  the 
hydrogen  entry  (charging,  x  =  0)  and  exit  sides  of  the  iron  (steel, 
x  =  L)  membrane.  C  and  C  are  the  hydrogen  concentrations  in  the 
membrane  atx=0  and  x=  L,  respectively. 


toward  the  paJladlum-coated  surface  where  it  gets 
oxidized  at  a  rate  that  is  measured  by  the  permeation 
current,  ip.  For  the  present  system,  it  has  been 
shown  that  the  rate  of  the  permeation  process  is 
controlled  by  diffusion  of  atomic  hydrogen  within  the 
membrane,  such  that: 


where  F  is  the  Faraday  constant,  D  is  the  dlffuslvlty 
of  hydrogen  within  the  membrane.  L  is  its  thickness, 
and  1  is  the  steady-state  value  of  ip. 

Figure  2  shows  the  effect  of  thiosulfate  concen-  . 
tration  on  the  hydrogen  permeation  transients 
through  the  steel  membranes  for  a  cathodic  charging 
rate  of  1.25  mA/cm^  in  the  0.05-M  H2SO4  +  0.45-M 
NajSO.  (blank)  solution.  Note  the  profound  effect  of 
the  concentration  of  thiosulfate  on  ip,  particularly  the 
1..  A  0.1 -mM  thiosulfate  addition  to  the  blank  solu¬ 
tion  Increases  the  L  by  about  twofold,  while  the  addi¬ 
tion  of  10  mM  thiosulfate  Increases  i_  by  about 
sixfold.  This  demonstrates  that  the  addition  of  thio¬ 
sulfate  Increases  the  1.  of  hydrogen  within  the  steel 
In  acidic  (pH  1.8)  solution.  This  is  attributed  to  an 
Increase  in  C".  which  results  from  the  effects  of  thio¬ 
sulfate  (and/or  its  decomposition  products)  on  the 
rate  of  hydrogen  absorption  within  the  membrane 
(Reaction  [3]).  Ixr'additlon,  thiosulfate  enhances  the 
HER  in  the  pH  1.8  solution  (Reactions  (1)  and  [2]). 
Figure  3  shows  Tafel  plots  of  the  HER  obtained  on 
steel  at  different  thiosulfate  concentrations  (-r)  = 

E  -  E"!.  E"i  =  -0.349  V  vs  saturated  calomel  electrode 
(SCEl).  This  figure  shows  that  an  increase  In  thiosul¬ 
fate  concentration  enhances  the  HER  (l.e.,  the  hydro¬ 
gen  oveipotential  decreases).  The  Tafel  slope 
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FIGURE  2.  Effect  of  SzOl~  concentration  on  the  hydrogen  permeation 
transients  obtained  on  a  steel  membrane  of  thickness  0.25  mm  in 
0.05  M  H2SO4  -h  0.45  M  Na2S04  (pH  1.8)  at  a  cathodic  charging 
current  density  of  1.25  mA/crrf. 


0.1  1  10 


ic  (mA/cm^) 

FIGURE  3.  Effect  of  concentration  on  the  Tafei  plots  of  the 
HER  obtained  on  a  steel  membrane  of  thickness  0.25  mm  in  0.05  M 
H2S04  +  0.45  M  Na2S04. 


Time  (s) 

FIGURE  4.  Effect  of  SzO^r  concentration  on  the  hydrogen  permeation 
transients  obtained  on  a  steel  membrane  of  thickness  0.25  mm  In 
0.5  M  NazSOM  (pH  7)  at  a  charging  current  density  of  1.25  mA/crrf. 


amounts  to  130  mV  =  2.3RT/0.45F  for  the  conditions 
in  Figure  3.  This  compares  well  with  a  theoretical 
value  of  2.3RT/0.5F  for  the  coupled  discharge- 
recombination  mechanism  which  is  to  be  expected 
when  the  Langmuir  isotherm  is  applicable. 


In  acidic  media,  thiosulfate  readily  decomposes 
to  give  colloidal  sulfur  and  sulfurous  acid  (H2SO3): 

S2OI''  +  2FF  S  +  H2SO3  (5) 

Consequently,  the  increase  in  the  rates  of  the  HER 
(Reactions  [1]  and  [2])  and  of  the  LIAR  (Reaction  [3]) 
may  be  attributed  to  the  remaining  thiosulfate  ions, 
colloidal  sulfur,  H2SO3,  and/or  the  reduction  prod¬ 
ucts  of  any  of  the  above  species. 

The  next  set  of  experiments  were  carried  out  on 
thiosulfate  in  a  neutral  medium,  where  Reaction  (5) 
is  not  likely  to  occur  and  the  predominant  species  in 
solution  is  the  thiosulfate  ion.  The  hydrogen  perme¬ 
ation  transients  in  the  neutral  medium  are  shown  in 
Figure  4.  This  figure  shows  that  thiosulfate  increases 
the  permeation  rate  of  hydrogen  within  iron  in  the 
neutral  medium.  A  concentration  of  0. 1  mM  thiosul¬ 
fate  increases  L  by  about  twofold  (which  is  compa¬ 
rable  to  that  observed  in  the  acid  medium)  while  the 
10-mM  concentration  increases  L,  by  about  threefold 
(which  is  less  than  that  observed  in  the  acid  me¬ 
dium).  This  increase  in  the  hydrogen  permeation  rate 
in  the  neutral  solution  was  also  accompanied  by  an 
increase  in  the  rate  of  hydrogen  evolution.  Figure  5 
shows  the  current-overpotential  plots  obtained  in 
0.5  M  Na2S04  at  different  thiosulfate  concentrations 
(t|  =  E  ~  =  -0.661  Vsce).  Hence,  thiosulfate  also 

depolarizes  the  HER  in  the  neutral  medium.  The 
Tafei  slope  is  200  mV,  2.3  RT/0.29F.  in  the  blank 
electrolyte  and  in  the  presence  of  0. 1  mM  and  1  mM 
thiosulfate.  This  indicates  that  the  mechanism  of  the 
HER  is  no  longer  coupled  discharge-recombination 
that  requires  a  Tafei  slope  of  120  mV  =  2.3  RT/0.50F. 
Note  tliat  in  the  10-mM  thiosulfate  solution  the  Tafei 
slope  is  120  mV.  However,  this  value  is  affected  con¬ 
siderably  by  the  cathodic  reduction  of  thiosulfate  at 
the  relatively  high  concentration  of  10  mM. 

Comparison  of  Figures  2  and  4  reveals  that  the 
rate  of  hydrogen  permeation  within  iron  is  greater  in 
an  acidic  than  neutral  medium,  for  both  the  blank 
solution  and  in  the  presence  of  a  given  amount  of 
thiosulfate  (or  its  decomposition  products).  The  in¬ 
crease  in  L,  produced  by  the  thiosulfate  is  also 
greater  in  the  acid  solution.  This  points  to  a  possible 
enhancing  effect  also  by  the  decomposition  products 
of  thiosulfate  on  hydrogen  permeation.  To  test  this 
hypothesis,  experiments  were  run  in  the  presence  of 
sodium  sulfite  (NasSOg)  in  the  acidic  medium  (pH  1.8) 
at  the  same  concentrations  as  in  the  thiosulfate  ex¬ 
periments  (Figure  6).  It  was  revealed  that  sulfite 
causes  an  increase  in  the  hydrogen  permeation  rate. 

This  is  also  coupled  by  an  increase  in  the  rate  of  the 
HER  (Figure  7).  As  a  result  of  the  proteolytic  equilib¬ 
ria  involving  sulfite  ions,  which  are: 

H2SO3 H^  +  HSO3  kj  =1.2x10“^  (6) 
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FIGURE  5.  Effect  of  S20t  concentration  on  the  Tafet  plots  of  the 
HER  obtained  on  a  steel  membrane  of  thickness  0.25  mm  in  0.5  M 
Na^SO,. 


FIGURE  6,  Effect  ofSO§~  concentration  on  the  hydrogen  permeation 
transients  obtained  on  a  steel  membrane  of  thickness  0.25  mm  in 
0,05  M  H2SO4  +  0.45  M  Na2S04  (pH  1.6)  at  a  charging  current  density 
of  1.25  mA/cm^. 


HSO3  +  SO^"  kj  =  6,2  X 10^^  (7) 

the  test  solutions  in  the  acid  medium  contain  pre¬ 
dominantly  H2SO3,  HSO5,  and  to  a  much  lesser  ex¬ 
tent,  SO§“.  In  the  neutral  medium,  where  the 
predominant  species  is  SOf'.  NasSOg  had  no  effect  on 
either  the  HER  or  the  hydrogen  permeation  current 
within  the  iron. 

The  extent  of  the  enhancement  in  the  hydrogen 
permeation  rate  at  different  thiosulfate  or  sulfite  con¬ 
centrations  may  be  expressed  quantitatively  by  an 
enhancement  factor,  e: 

i«  (additive) 

where  L  (additive}  and  L  (blank)  are  the  steady- state 
hydrogen  permeation  current  densities  for  the  solu¬ 
tion  with  and  without  the  additive,  respectively. 

Table  1  lists  values  of  L  and  e  for  different  concentra¬ 
tions  of  thiosulfate  and  NaaSOs  in  acidic  and  neutral 
solutions.  The  table  reveals  that  thiosulfate  (and/or 
its  decomposition  products)  In  acid  medium  shows  a 
higher  promoting  effect  (larger  e)  than  NaaSOg  at  the 
same  concentration.  If  one  assumes  that  1  mol  of 
(SOf"  +  HSO3’  +  H2SO3)  results  from  the  decomposi¬ 
tion  of  1  mol  of  thiosulfate,  it  can  be  reasoned  that 
the  resulting  colloidal  sulfur  and/or  its  reduction 
products  also  cause  some  enhancement  of  the  hydro¬ 
gen  permeation.  In  addition,  thiosulfate  was  found 
to  promote  hydrogen  permeation  within  iron  in  the 
neutral  medium. 

In  addition  to  their  effects  on  hydrogen  evolution 
and  permeation  within  iron,  SaOi",  SOI",  and  S  are 
thermodynamically  unstable  at  the  rather  negative 


FIGURE  7.  Effect  of  SO|-  concentration  on  the  Tafel  plots  of 
the  HER  obtained  on  a  steel  membrane  of  thickness  0.25  mm 
in  0.05  M  HzSOi  +  0.45  M  Na2S04  at  a  charging  current  density  of 
1.25  mA/cnf. 


potentials  of  these  experiments^®  and,  therefore,  may 
be  reduced  on  the  metal  surface  to  hydrogen  sulfide 
(H2S).  H2S  is  well  known  to  enhance  hydrogen  ab¬ 
sorption  into  metals  and  to  enhance  the  HER.^^-^® 
However,  as  is  shown  in  the  literature,  these  reduc¬ 
tion  reactions  do  not  play  a  significant  role  under  the 
conditions  of  these  experiments,®^  Even  very  small 
amounts  of  H2S.  however,  could  contribute  to  the 
observed  enhanced  hydrogen  permeability.®® 
y 

CONCLUSIONS 


❖  Thiosulfate  promotes  hydrogen  evolution  and  per¬ 
meation  within  iron  In  acidic  and  neutral  solutions. 
Its  effects  in  the  acid  medium  are  attributed  mainly 
to  Its  decomposition  products  (I.e.,  H2SO3,  HSOi.  SOi, 
and  colloidal  sulfur)  since  thiosulfate  is  readily  de- 
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TABLE  1 

Effects  of  the  Concentrations  ofS^Ot  andSOt  on  the  Hydrogen  Perrneation  Rate  within  Iron 
and  the  Enhancement  Factor  (e)  in  Acidic  (pH  1.8)  and  Neutral  (pH  7)  Solutions  at  a  Charging  Current  of  1.25  mA/cnf 


c 

(mM) 


_ pH1»8 _ _ 


SOI' 

_ pH  1.8 _ 

L  (uA/cm^)  e 


s^oi- 

_ pH^ _ 

I*  (pA/cm®)  E 


0 

0.1 

1 

10 


3.2 

1 

3.2 

6.6 

2.0 

5.9 

14.0 

4.3 

12.4 

20.4 

6.3 

16.8 

1 

.1.8 

3.8 

5.2 


0.8 

1 

1.3 

1.6 

1.7 

2.1 

2.2 

2.8 

composed  in  acidic  media  to  give  the  above  products, 
and  since  results  with  only  sulfite  in  the  same  acidic 
solution  show  sHghtly  less  enhancement  of  hydrogen 
permeation.  Thiosulfate  in  neutral  solution,  where  it 
remains  mainly  undecomposed,  also  causes  en¬ 
hancement  of  the  hydrogen  absorption  reaction  but 
to  a  lesser  extent.  On  the  other  hand,  sulfite  in  the  . 
neutral  medium  has  no  effect  on  either  process. 

Thus,  the  major  effects  of  thiosulfate  are  caused  by 
Its  decomposition  product  in  acidic  media  and  to 
thiosulfate,  itself,  in  neutral  media.  The  mechanisms 
involved  in  the  enhancement  processes  are  yet  to  be 
determined. 
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